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ABSTRACT 

A study is made of a sample of 58 dust-enshrouded Asymptotic Giant Branch 
Q>^ i (AGB) stars (including 2 possible post AGE stars), of which 27 are carbon-rich 

and 31 are oxygen-rich. These objects were originally identified by Jura & Klcin- 
mann as nearby (within about 1 kpc of the sun) AGB stars with high mass-loss rates 
(M > 10~^ M0?/r~^). Ground-based near-infrared photometry, data obtained by the 
IRAS satellite and kinematic data (radial and outflow velocities) from the literature 
are combined to investigate the properties of these stars. The light amplitude in the 
near-infrared is found to be correlated with period, and this amplitude decreases with 
increasing wavelength. Statistical tests show that there is no reason to suspect any 
' difference in the period distributions of the carbon- and oxygen-rich stars for periods 

less than 1000 days. There are no carbon-rich stars with periods longer than 1000 days 
in the sample. The colours are consistent with those of cool stars with evolved circum- 
stellar dust-shclls. Luminosities and distances are estimated using a period-luminosity 
relation. Mass-loss rates, estimated from the 60 fxm fluxes, show a correlation with 
various infrared colours and pulsation period. The mass-loss rate is tightly correlated 
O ' with the K — [12] colour. The kinematics and scale-height of the sample shows that 

the sources with periods less than 1000 days must have low mass main-sequence pro- 
^ ■ genitors. It is argued that the three oxygen-rich stars with periods over 1000 days 

^ ' probably had intermediate mass main-sequence progenitors. For the other stars an 

average progenitor mass of about 1.3 is estimated with a final white dwarf mass 
of 0.6 Mq. The average lifetime of stars in this high mass-loss AGB phase is estimated 
rS I to be about 4 x 10^ years, which suggests that these stars will undergo at most one 

■ more thermal pulse before leaving the AGB if current theoretical relations between 

thermal interpulse-period and core-mass are correct. 
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1 INTRODUCTION 

Stars of low to intermediate main-sequence mass, up to 
about 6-8 Mq , will eventually evolve up the Asymptotic Gi- 
ant Branch (AGB) (Iben & Renzini 1983). At the tip of 
the AGB, where they combine a low effective temperatures 
(Te// ~ 3000 K) with high luminosities (L ^ 3000 L©), the 
stars pulsate and are subject to rapid mass loss. Typical ex- 
amples of such stars are the Miras and the oh/ir[] stars 
(e.g. Feast & Whitelock 1987; Habing 1996). Their light 
output varies on timescales from around one hundred days 

* A small fraction of OH/IR stars are supergiants, i.e. their pro- 
genitors were massive stars. 



to well over a thousand days. They also lose mass rapidly 
(10"'' ^ M $ 10'* MQyr'^). They follow various period- 
luminosity relations and are at or near the maximum lu- 
minosity they will ever achieve during their evolution. This 
makes them very useful as distance indicators and as probes 
of galactic structure. They are also important contributors 
of dust into the interstellar medium, which may be crucial 
for the formation of new stars and planets. 

Jura & Kleinmann (1989 henceforth JK89) identified a 
group of mass losing AGB stars within about 1 kpc of the 
sun. In 1989 a monitoring programme of stars in the JK89 
sample visible from the southern hemisphere was initiated 
from SAAO. This entailed obtaining near-infrared photo- 
metry at JHKL over a time interval of several years, from 
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which pulsation periods could be determined. The aim of 
this study was to investigate the physical properties of the 
sample, in particular variability, colours, mass-loss rates and 
kinematics, using SAAO photometry together with infrared 
photometry and kinematic data from the literature. 

Finally these data are used to estimate the average pro- 
genitor mass and lifetime of stars in this high mass-loss 
phase. 



2 THE SAMPLE 

JK89 aimed to develop a more quantitative understanding of 
mass loss from local AGB stars. They identified AGB stars 
undergoing heavy meiss loss (M > 10~^ Mq j/r"^), hence- 
forth referred to as dust-enshrouded, within about 1 kpc of 
the sun, and within -33° < 5(1950) < 82°. Their selection 
was based on infrared photometry from the Two Micron Sky 
Survey, IRAS and the AFGL survey from which they iden- 
tify 63 sources. They give reasons for thinking this sample 
was reasonably complete. 

A few comments should be made at this point. JK89 
assume a luminosity of 10* Lq for all of these AGB stars. 
It will be seen later that, if the period luminosity relation 
(PL) is correct, this leads to the inclusion of more luminous 
OH/IR stars outside the 1 kpc sphere and up to ~ 2 kpc 
away. They note that some red supcrgiants will satisfy their 
selection criteria; those they largely excluded, but included 
the star 3 Puppis, an A supergiant, with the idea that it 
might be undergoing the transition from AGB star to plan- 
etary nebula. However, Plets, Waelkcns & Trams (1995) con- 
vincingly showed that 3 Puppis is a genuine massive super- 
giant, and we therefore exclude it from this study. JK89 also 
included the carbon-rich planetary nebulae NGC 7027 (e.g. 
Justtanont et al. 2000) in their sample; this object is ex- 
cluded from the present study. The remaining 61 sources are 
listed in Table 1 (this includes 3 M supergiants, 17513 — 2313, 
18135 - 1641 and 18204 - 1344 (Ukita & Goldsmith 1984; 
Winfrey et al. 1994) of which we only became aware part the 
way through this project) . Near-infrared data are tabulated 
for two of them, but none are used in the analysis. 

Chemical types (O and C for oxygen- and carbon-rich, 
respectively) for all the sources were taken from Loup et 
al. (1993), and spectral types obtained from the SIMBAD 
data base for 51 of the sources are also listed in Table 1. 
Although 04307 + 6210 is listed as M6 in SIMBAD, Groc- 
newegen (1994) suggests this is incorrect and the star is 
actually carbon rich; therefore no spectral type is assigned 
to it here. The supergiants are marked (I). 

The final sample thus consists of 27 carbon-rich stars 
and 31 oxygen-rich stars (excluding the 3 supergiants) of 
which 2 are S-type. Almost all of the stars are of very late 
spectral type (M, C, N, R and S), with two exceptions: the 
Red Rectangle and the Egg Nebula. Both are well-known 
carbon-rich bipolar reflection nebulae (see, e.g. Sahai et al. 
1998; Waters et al. 1998). The Red Rectangle is a wide bi- 
nary system with am oxygen-rich dust disk (Waters et al. 
1998). These two objects are generally thought to be post- 
AGB stars on their way to become Planetary Nebulae. 

The Benson et al. (1990) catalogue lists sources which 
have been examined for OH, SiO and H2O maser emission. 
Of the 28 oxygen-rich stars in our sample which have been 



Table 1. Infrared Sources 



INFRARED 


CHEM. 


SPEC. 


VARIABLE 


NO. OI 


SOURCE 


TYPE 


TYPE 


NAME^ 


REF.* 


IRAS ()()()42 + 4248 
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MIO 


KU And 


47 
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M9 
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216 


TR A*^ ni 1 eiQ _|_ 799n 
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Se 


S Cas 


57 


IRAS 02270 — 2619 






R For 


91 


TR A*^ n9^1 ft -1- PiA'^'^ 
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M9 


V UUU V^OjO 


39 


IRAS 02351 — 2711 


o 


M9 


UU For 


37 


IRAS 03229 + 4721 






V384 Per 


101 


TR A ^ O'^'^n? -1- 1 1 1 ""i 


Q 


M6e 


IK Tau 


277 


IRAS (J43U7 + 6210 








38 


IRAS (J4566 + 5606 


Q 


MS. 5 


TX Cam 


147 




o 


M 


NV Aur 


120 


IRAS 05411 + 6957 


o 


M9 


BX Cam 


52 


TR A*^ n^^^Q -1- 74^n 


o 


M7 


V Cam 


72 


TR AS nfil 7R 1 n^fi 


Q 


B8V 




272 


IRAS 06300 + 6058 


o 


IVI9 


AP Lyn 


74 


IRAS 06500 + 0829 


Q 


M9 


GX Mon 


52 


IRAS 08088 — 3243 


Q 




V346 Pup 


31 


IRAS 09116 — 2439 








27 


IRAS 09429 — 2148 


Q 


M9 


IW Hya 


74 


TR A 1^ HQd'^O -i- 1 '^'^n 


Q 


Q 


VV J-JCU 


1004 


TR AS 1 ni ^1 -I- ^fl/lQ 


Q 


Ce 


RW LMi 


218 


TRAS 1 n4.Qi 9n'=;Q 


Q 


Q 


V iiya 


164 


IRAS 12447 + 0425 


Q 


Re 


RU Vir 


82 


IRAS 17049 — 2440 


Q 






37 


IRAS 17119 + 0859 


Q 


M9 


V2108 Oph 


51 


IRAS 17297 + 1747 


o 


M2 


V833 Her 


50 


IRAS 17360 — 3012 


Q 




V1019 Sco 


40 


IRAS 17411 — 3154 


Q 






31 


TR AS 17'^1'^ — 9*^1^ 


o 


M5(I) 


V774 So-r 


15 


TRAS 1 SnnQ 901 Q 


o 


MS 


V/11 90 So-r 


43 


TRAS 1 sn4n — nQ4i 


Q 


Q 


FX Ser 


43 


IRAS 18135 1641 


o 


M5(I) 




15 


IRAS 18194 — 2708 








33 


IRAS 18204 — 1344 


Q 


M8(I) 




39 


IRAS 18240 + 2326 


c 


Q 




51 


IRAS 18333 + 0533 


() 




NX Scr 


44 


IRAS 18348 — ()r)2f) 


() 


M 


V437 Set 


164 


TR A S 1 8'^zLQ -1- 1 09'-! 


Q 


M9 


VI 111 OnH 

V J- J- J- J- \j \j 11 


85 


TR A S 1 S^Q7 -1- 1 


Q 


Q 


V821 Her 


75 


IRAS 18398 — 0220 






V1417 Aql 


46 


IRAS 18413 + 1354 


o 


M7 


V837 Her 


35 


TR A S 1 R^fin — 9Q^d. 

1 1 V -r\ O XOUL*L/ — ik?U^ 


o 


M9 


V^Q'i^ So-r 
V Ot/OO 


28 


IRAS 19008 + 0726 




R 


V1418 Aql 


43 


IRAS 19059 — 2219 


Q 


M8 


V3880 Sgr 


53 


IRAS 19093 — 3256 


Q 


M9 


V342 Sgr 


20 


IRAS 19126 — 0708 


Q 


S 


W Aql 


97 


TR AS 1 Q1 7'^ — 0807 


Q 


Q 


V1420 Aql 


52 


TRAS ^Q^^9^ -i- 97^17 


Q 


Q 




48 


TRAS 90077 — 0fi9^ 


o 


M 


VI ^nn A(-ii 


73 


IRAS 20396 + 4757 




N 


V Cvs 


162 


IRAS 20440 - 0105 


o 


M9 


FP Aqr 


30 


IRAS 20570 + 2714 


c 


Ce 




35 


IRAS 21032 - 0024 


c 


Cc 


RV Aqr 


54 


IRAS 21286 + 1055 


o 


M7c 


UU Peg 


35 


IRAS 21320 + 3850 




C 


V1426 Cyg 


72 


IRAS 21456 + 6422 


o 


M6 


RT Cep 


18 


IRAS 23166 + 1655 


c 




LL Peg 


109 


IRAS 23320 + 4316 


c 


C 


LP And 


98 


IRAS 23496 + 6131 


o 


M9 


V657 Cas 


24 


RAFGL 1406 


c 




lY Hya 


40 


RAFGL 2688 


c 


F5Iae 


(Egg Nebula) 


347 



f-Namcs in brackets are most commonly used. 
*-From 1960 to 1999. 

examined for OH and H2O, 20 were detected in OH and 20 
in H2O. While 19 of the 21 sources searched for SiO were 
detected. Interestingly, 4 of the 11 carbon-rich stars which 
were searched for SiO also had positive detections. We can 
thus estimate that 70 percent of the oxygen-rich stars have 
OH and H2O maser emission and about 90 percent have SiO 
emission. 
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Table 2. Near-infrared data. The Date is given as JD-2440000.00. 



DATE 


J 


H 


JV 


T 

±J 
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J 


H 


K 


L 


IRAS 01037 + 1219 






IRAS 01037+ 1219 .... 


.continued 


06304 57 


7.04 


3.90 


1 70 

-L. 1 U 


— u.oo 


08589 


8.71 


5.32 


3.04 


0.58 


n63'^4 4Q 


6.98 


3.84 


-L .Uo 


— u.uy 


08855 58 


6.60 


3.75 


1.74 


—0.58 




7.00 


3.84 




— u.oy 


UOO ( '±,00 


6.57 


3.70 


1.71 


—0.59 




6.99 


3.86 


1 79 


n fi7 


080*^0 A'^ 


6.51 


3.64 


1.67 


—0.60 


L/UOf?^ .OU 


7.08 


3.92 


1 7fi 


n fizi 

— U.D^ 


uoyu J- . oi 


6.52 


3.63 


1.68 


—0.56 




7.27 


4.05 


1 Q9 


n 

— U.^D 


O C/ Cj O . ij 


6.57 


3.70 


1.76 


—0.52 


06640 67 


9.66 


5.84 




U.DO 


0921 2 61 


8.40 


5.22 


2.99 


0.50 


06657 61 


9.81 


5.90 


o.oU 


U.Do 


09232 54 


8.52 


5.30 


3.03 


0.55 


06694 55 


9.80 


5.79 






09300 40 


8.68 


5.43 


3.06 


0.41 


0671 3 45 


9.79 


5.93 




u.oo 


09582 57 


6.19 


3.34 


1.36 


—0.81 


06741.41 


9.74 


5.89 


^ 9^ 


A7 


09614.53 


6.30 


3.43 


1.45 


—0.73 


06750 38 


9.76 


5.90 


9fi 


n f^9 


09637.43 


6.41 


3.53 


1.54 


—0.62 


06754 33 


9.71 


5.90 


9R 


n 47 


09672.36 


6.65 


3.73 


1.72 


—0.42 


06774 33 


9.57 


5.78 


O. -LO 


u.oy 


10004.45 


8.55 


5.11 


2.65 


0.15 


06775 31 


9.81 


5.75 


0. 1 i 


u.oo 


10019.44 


7.87 


4.67 


2.33 


—0.10 


06784 28 


9.23 


5.52 


9 Q"^ 


n 9*^ 


10063.33 


6.74 


3.85 


1.75 


—0.55 


06801 27 


8.28 


4.88 


2 49 


— u.uy 


10089.28 


6.66 


3.76 


1.64 


—0.58 


07005 66 


6.78 


3.72 


-L.UO 


— u.oy 


10298.63 


6.37 


3.42 


1.44 


—0.59 


07024 55 


6.82 


3.74 


1.69 


—0.60 


10360.49 


6.88 


3.83 


1.79 


—0.23 


07072 4Q 


7.00 


3.90 


1 84 


n 47 

— u.^ / 


10417.35 


7.63 


4.38 


2.20 


0.14 


07114.40 


7.32 


4.17 


2 12 


—0 22 


10717.50 


7.03 


3.89 


1.70 


—0.56 


071 27 36 


7.46 


4.29 


2 25 


n 07 

— U.U ( 


10736.46 


6.83 


3.81 


1.63 


—0.62 


071 A7 0*^ 
U ( 1 .Uo 


i .Do 




9 '^f^ 
Z.oO 


n c\A 
u.u^ 


10805.30 


6.39 


3.37 


l.OU 


n SI 


07367.67 


9.94 


6.07 


'^A 


u.ou 












07461.42 


8.38 


5.02 


9 i^l 
Z.Di 


n m 
— u.u± 


IRAS 02270 - 2619 






07735.60 


7.03 


3.97 


1.98 


—0.35 


10361.45 


4.64 


2.88 


1.54 


0.21 


07792.51 


7.60 


4.40 


2.35 


0.01 


10399.48 


4.49 


2.87 


1.45 


0.07 


07817.42 


7.93 


4.63 


2.50 


0.14 


10437.32 


4.20 


2.54 


1.27 


-0.10 


07835.38 


8.15 


4.82 


2.64 


0.25 


10467.28 


3.85 


2.22 


1.05 


-0.30 


08109.61 


8.21 


4.80 


2.40 


-0.21 


10499.25 


3.70 


2.06 


0.92 


-0.38 


08164.51 


6.87 


3.94 


1.87 


-0.53 


10673.63 


4.72 


2.90 


1.54 


0.12 


08211.34 


6.66 


3.74 


1.72 


-0.60 


10754.38 


4.40 


2.64 


1.34 


0.08 


08232.37 


6.53 


3.67 


1.67 


-0.61 


10796.35 


4.18 


2.51 


1.25 


-0.03 


08473.66 


7.59 


4.45 


2.42 


0.18 


10854.29 


3.49 


1.94 


0.83 


-0.37 


08495.59 


7.84 


4.63 


2.57 


0.30 












08520.53 


8.09 


4.87 


2.73 


0.41 


IRAS 02351 - 2711 






08554.43 


8.46 


5.06 


2.89 


0.36 


09170.68 


3.32 


1.99 


1.24 


0.44 



DATE J H K L 



DATE J H K L 
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u. 
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0.05 
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\ 
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0.08 
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1 
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7S 
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I 
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0.27 
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1 
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2 
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0.40 
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1 
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os8fin sn 


2, 


.18 


0.59 
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]^ 


A^ 
.^0 


UOOUO-O 1 


2 
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0.54 


n 4n 


1 

— J., 


t;i 
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nS868 S7 


2 


.25 


0.61 


— u.oo 


_][ 


47 


03879 35 


2 


.39 


0.72 


9^^ 


1 

— 1 , 


40 
.4U 


KjOO<y^ . OZ 


2 


.55 


0.81 


91 


1 

— 1 , 


.00 


uooyu . 00 


2 


.51 


0.82 


1 Q 

— u. ±y 


1 

— 1 , 


.oy 


uiuyu . 00 


\ 


.13 


—0.09 


QA 


9 


1 7 


04101 iS9 




.09 


—0.18 


1 00 
1 .uu 




1 Q 

. ly 


041 58 58 




.05 


—0.32 


110 
— 1 . lU 


_2 


1 Q 

. ly 


041 fil ^9 
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—0.32 


1 in 


9 


99 

.zz 


(1494^ QO 


1 

J., 


.00 


n no 


— U.OCJ 


^ 


87 
.0 J 


04252.31 


1, 


.59 


0.04 


n SI 


1 

— 1 , 


.00 


04284.28 


1, 


.90 


0.28 


fii 


_^ 


.uo 


04578.42 


1, 


.32 


-0.15 


1 

— i .uo 


9 

— z, 


9^^ 
.ZO 


04582.36 


1, 


.34 


-0.09 




— 2, 


.24 


04604.36 


1, 


.32 


-0.19 


-1.08 


-2, 


.15 


04605.34 


1, 


.25 


-0.22 


-1.12 


-2, 


.32 


04607.32 


1, 


.26 


-0.24 


-1.18 


-2, 


.36 


04614.33 


1, 


.27 






-2, 


.32 


04623.33 


1, 


.27 






-2, 


.29 


04627.30 


1, 


.37 






-2 


.23 


04634.29 


1 


.31 






-2, 


.25 


04911.46 


2, 


.81 


0.89 


-0.32 


-1, 


.66 


04914.43 


2, 


.68 


0.80 


-0.41 


-1, 


.81 


04918.45 


2, 


.73 


0.85 


-0.34 


-1, 


.73 


04947.38 


2, 


.50 


0.77 


-0.45 


-1, 


.76 


04949.40 


2, 


.48 


0.72 


-0.43 


-1, 


.72 


04954.40 


2, 


.32 


0.68 


-0.48 


-1, 


.78 


04958.34 


2, 


.19 


0.55 


-0.57 


-1, 


.84 
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04960 .S5 


2.15 


0.56 


—0.56 


-1.84 


0521 8 61 


1.87 


0.18 


—0.71 


1 77 


yJO^O 1. . OU 


2.24 


0.49 


—0 52 


1 fin 


yjoooo 'Oo 


2.73 


0.87 


—0 30 


-1.56 


n'^(^09 fi9 
uouuz. uz 


1.25 


—0.30 


— 1.14 


9 9c; 
— Z.ZD 


n^^fil 9 ^^7 

UOU LZi.O 1 


1.33 


—0.26 


— 1.11 


9 97 
— Z.Z ( 




1.33 


—0.23 


— 1.09 


9 1 fl 
— Z. lU 


05647 52 


1.58 


—0.05 


—0.94 


9 09 
— Z.UZ 


05681 34 


1.92 


0.20 


—0.75 


— 1.00 


0571 30 


2.29 


0.50 


—0.51 


-1.58 (-~| 

— i.OO 

— 1.59 ^ 


06309.58 


2.81 


0.89 


—0.28 


06335 59 


2.68 


0.83 


—0.32 


06378 51 


1.75 


0.24 




1 QQ CT) 

— ±.yy 


06439 29 


1.27 


0.00 




2 23 

9 "^Ifi 
z.oo ^ 


06491 25 


1.16 


—0.40 


— 1.26 


06645 64 


2.46 


0.54 


—0.52 


-1.62 2 


06658 65 


2.65 


0.71 


—0.40 


— 1 55 S-^ 
en 

-1.41 


06690 59 


3.07 


1.02 


—0.17 


06723 57 

\J\J 1 £iKJ' kJ I 


3.28 


1.13 


—0.13 




06752 44 


3.06 


0.98 


—0.23 


— 1.01 


06777 37 


2.73 


0.79 


—0.38 


1 7n to 


06783 37 


2.73 


0.76 


—0.41 


1 fi7 ^ 


06803.30 


2.-58 


0.68 


-0.45 


-1.76 ^ 
-1.81 ^ 


06826.28 


2.38 


0.56 


—0.52 


07014.69 


l.dO 


-0.24 


-1.13 


— 2 20 


07056.60 


1.77 


0.04 


-0.89 


-1.97 r+ 


07074.61 


2.00 


0.24 


-0.72 


-1.90 ^ 


07114.54 


2.60 


0.68 


-0.38 


-1.56 Co 


07144.36 


2.96 


1.00 


-0.14 


-1.42 ^ 


07186.32 


3.23 


1.17 


-0.05 


-1.40 ^ 


07386.64 


1.42 


-0.10 


-1.04 


-2.26 1^ 
-1.93 ^ 


07506.38 


1.88 


0.14 


-0.83 


07530.34 


2.00 


0.26 


-0.73 


-1.87^' 


07805.58 


1.55 


-0.01 


-0.99 


-2.23 §1 


07821.50 


1.45 


-0.10 


-1.05 


-2.30 1 


07841.45 


1.32 


-0.22 


-1.15 


-2.37 3. 

-2.46 



07868.45 


1.21 


-0.36 


-1.27 



CO 



Table 2. ...continued. Near-infrared data. The Date is given as JD-2440000.00. 



DATE 


J 


H 


K 


T 
i J 


DATE 


J 


H 


i\ 


T 


IRAS 03507 + 1115 






IRAS 06176 - 1036 






07905 32 


1.21 


-0.42 


—1.33 


2 46 


06867.33 


6.61 


4.95 


^0 


l.OU 


07010 '?2 


1.23 - 


-0.42 


— 1.32 


—2 45 


10437.52 


6.70 


5.03 


3 44 




n^^l 9Q fi9 


3.42 


1.24 


—0.06 


1 A"^ 
— i .^o 












UO lUO .oo 


3.25 


1.11 


—0.17 


1 f^i 


IRAS 08088 - 3243 






n^^99S 41 

y}<j^Zj(j .'-I J. 


3.15 


1.10 


—0.16 


1 f^Q 

— i . oy 


09117.32 


9.67 


6.64 


A 9A 


1 AS 


nS2^2 41 


2.52 


0.70 


—0.47 


1 

— 1 .o^ 


09734.55 


8.95 


6.35 


A 90 

^. zu 


i-.OO 


L/O O L/ . O iJ 


1.63 


0.09 


—0.87 


9 1 A 


09832.29 


8.77 


6.19 


A HQ 


1-OD 


0S510 fiO 

VJOiJ i-tJ .\J\J 


2.76 


0.78 


—0.33 


1 AS 


10034.53 


7.43 


5.01 


09 


U.Oo 


08630 34 


4.09 


1.80 


0.39 


]^ ^2 


10090.52 


7.40 


4.98 


'i ni 

O.U J. 


U.tJ 1 


08873 fi4 


1.87 


0.16 


—0.87 


9 1 


10112.44 


7.46 


5.04 


o.uo 




08933 48 


2.37 


0.54 


—0.53 


1 70 


10478.46 


8.62 


6.10 




1 40 


08961 40 


2.75 


0.84 


—0.29 


— 1 52 


10500.38 


8.42 


5.92 


o.oy 


1 '^n 


089Q0 38 


3.19 


1.21 


—0.01 


1 'iA 

— i .o^ 


10591.24 


7.45 


5.07 


1 9 

o. iz 


U.DD 


0921 2 fiS 


2.59 


0.82 


—0.41 


— 1 92 


10796.45 


8.85 


6.18 


A 


1 40 


09236 62 


2.25 


0.56 


—0.61 


9 n*^ 


10916.39 


9.73 


6.97 


A 71 


1 Q1 


09581 67 


3.88 


1.64 


0.21 


1 '^l 


10978.21 


9.20 


6.51 


^.O J. 


1 62 


09614 60 


3.82 


1.64 


0.21 














09637 56 


3.89 


1.63 


0.21 


— 1 "^4 


IRAS 09116 - 2439 






09672 42 


2.80 


1.03 


—0.24 


— 1 7f\ 


07175.50 




10.93 


7.36 


3.03 


09709 34 


2.41 


0.66 


—0.55 


9 


09117.34 


13.93 


10.31 


f\ 77 




uy 1 zo.oo 


2 34 


u.ou 


— u.uo 


— 2 10 


09734.59 




9.75 


u.oy 


2 40 


09986.61 


3.79 


1.49 


0.12 


— 1 20 


09829.27 




10.21 


SQ 
u.oy 


9 Q1 


10019.52 


3.89 


1.56 


0.14 


1 97 


09887.23 




10.39 


7 07 


o. i 1 


10052.47 


3.68 


1.41 


0.01 


— i.oy 


10034.56 




9.15 


O.oo 


i.yi 


10086.35 


3.52 


1.35 


-0.04 


-1.51 


10086.56 




8.64 


5.34 


1.51 


10109.29 


3.48 


1.33 


-0.05 


-1.54 


10111.46 




8.47 


5.18 


1.36 


10125.27 


3.09 


1.12 


-0.21 


-1.69 


10123.50 




8.43 


5.15 


1.33 


10362.54 


2.51 


0.50 


-0.61 


-1.71 


10222.31 




8.58 


5.27 


1.42 


10417.42 


3.12 


0.98 


-0.24 


-1.38 


10256.20 




8.74 


5.46 


1.65 


10437.38 


3.33 


1.16 


-0.10 


-1.33 


10438.58 




9.96 


6.68 


2.65 


10469.32 


3.52 


1.32 





-1.31 


10483.45 




10.26 


6.88 


2.80 


10503.00 


3.40 


1.23 


-0.11 


-1.44 


10590.25 




10.16 


6.86 


2.78 


10721.61 


1.80 


0.00 


-1.04 


-2.17 


10795.51 




8.62 


5.31 


1.46 


10753.45 


1.90 


0.04 


-1.00 


-2.12 


10913.38 


12.46 


8.65 


4.71 


1.92 


10805.38 


2.16 


0.26 


-0.82 


-1.91 


10978.23 




9.03 


5.77 


1.97 



DATE 


J 


H 


JV 


T 
±j 


DATE 


J 


H 


K 


L 


IRAS 09429 - 2148 






IRAS 09429-2148 


.continued 


06356 62 


6.16 


4.41 


OQ 
o.uy 


1 46 


08024.23 


4.79 


3.14 


1.80 


0.09 


06427 60 

\J\J'-t^ 1 .\J\J 


6.43 


4.55 


O. J-O 


1 46 




4.91 


3.28 


1.93 


0.18 




6.46 


4.59 




1 44 


0^^91 ^7 
uoz ±o . O i 


5.80 


4.04 


2.66 


1.01 




6.47 


4.58 




1 A'i 


uozou . oo 


6.07 


4.28 


2.85 


1.17 


nfi4K7 4*^ 


6.54 


4.64 


99 


1 4S 


0^^9^^1 '^O 

UOZO ± . Ok) 


6.22 


4.38 


2.94 


1.26 




6.39 


4.56 


0. 10 


1 "^S 
±.oo 




6.26 


4.44 


2.98 


1.27 


06^06 37 


6.35 


4.53 


•^19 
O. 


1 '^4 


08327 38 

V700iJ j .OO 


6.39 


4.52 


3.05 


1.31 


06541 36 


5.34 


3.80 


9 "^9 
Z.OZ 


74 


08377 26 


6.44 


4.57 


3.06 


1.31 


06568 21 


5.05 


3.59 


9 '^'^ 
z.oo 


n "^7 


08616 60 


4.79 


3.10 


1.75 


0.08 


06747 63 


4.78 


3.19 


1 Q7 


n '^4 


0861 7 59 


4.78 


3.09 


1.75 


0.11 


06778 52 


4.85 


3.27 


9 0^1 


41 


08639 52 


4.75 


3.05 


1.72 


0.08 


06808 54 


4.95 


3.37 


2 15 


4Q 


08665 47 


4.81 


3.08 


1.76 


0.17 


06826 47 


5.06 


3.47 


9 94 


U.Oo 


08701 38 

uo I \j i- . oo 


4.89 


3.19 


1.87 


0.27 


06834 43 


5.10 


3.50 


2 28 


fi9 
u.uz 


08721 29 

UO 1 ^ i- . ^ij 


5.01 


3.30 


1.96 


0.34 


OfiS47 47 


5.22 


3.60 


Z.oD 


7^ 
U. I 


08761 21 

UO 1 U J- . J- 


5.20 


3.46 


2.14 


0.56 


06873 39 


5.38 


3.76 


2 52 


S^i 
U.OO 


08992 53 


6.65 


4.54 


2.92 


1.17 


06894 33 


5.52 


3.91 


2 64 




u.yo 


09106.27 


5.68 


3.85 


2.36 


0.60 


06897 36 


5.49 


3.90 


9 fi^ 


u.yy 


09382.48 


5.17 


3.35 


2.01 


0.43 


06897 37 


5.57 


3.93 


2.68 


1.00 


09405.46 


5.31 


3.46 


2.10 


0.55 


06934 28 


5.82 


4.14 




1 1 

J.. lU 


09467.39 


5.78 


3.83 


2.41 


0.90 


06978 23 


6.17 


4.37 


09 

O.UZ 


-L.OU 


09501.25 


6.01 


4.01 


2.56 


0.98 


07121 55 


6.37 


4.50 


07 

O.U ( 


1 29 


09677.55 


6.73 


4.45 


2.79 


1.04 


071 71 


^ fiO 

o.ou 


'i QO 

o.yu 


Z.OO 


71 
U. 1 1 


09824.34 


4.81 


3.00 


i .DO 


OS 


07179.50 


5.37 


3.78 


2 44 


fi9 
U.UZ 


09887.27 


4.67 


2.82 


1.51 


-0.02 


07192.49 


5.21 


3.66 


Z.OO 


f^l 
U.Ol 


10034.58 


5.06 


3.17 


1.87 


0.43 


07213.49 


5.03 


3.51 


2.21 


0.41 


10088.59 


5.41 


3.47 


2.13 


0.71 


07236.38 


4.91 


3.38 


2.08 


0.29 


10463.56 


4.61 


2.83 


1.57 


0.06 


07268.32 


4.80 


3.25 


1.95 


0.20 


10500.44 


4.56 


2.75 


1.50 


0.00 


07533.54 


5.41 


3.76 


2.46 


0.84 


10800.47 


6.17 


4.06 


2.60 


1.05 


07542.41 


5.51 


3.82 


2.51 


0.93 












07581.42 


5.78 


4.04 


2.70 


1.07 


IRAS 09452 + 1330 






07617.33 


6.03 


4.24 


2.84 


1.20 


07215.45 


8.38 


5.01 


2.10 


-1.71 


07690.20 


6.56 


4.57 


3.05 




07242.36 


8.24 


4.85 


1.95 


-1.72 


07890.54 


5.06 


3.46 


2.10 


0.28 


07535.56 


6.39 


3.05 


0.28 


-3.23 


07895.59 


4.98 


3.41 


2.05 


0.20 


07583.42 


6.74 


3.40 


0.61 


-3.01 


07904.56 


4.95 


3.36 


2.01 


0.20 


07617.35 


7.05 


3.69 


0.90 


-2.67 


07921.48 


4.85 


3.27 


1.92 


0.11 


07669.24 


7.45 


4.12 


1.33 


-2.29 



Table 2. ...continued. Near-infrared data. The Date is given as JD-2440000.00. 



DATE 


J 


H 


K 


L 


DATE 


J 


H 


K 


T 
Ij 


IRAS 09452 + 1330 .... 


.continued 


IRAS 10491 - 2059 






07903 55 


7.70 


4.48 


1.70 


—1.94 


06486 41 


1.86 


0.32 


—0.72 


1 QQ 

— J..OO 




6.83 


3.52 


0.73 


—2.87 


nfi'=;n9 ao 


1.98 


0.42 


—0.64 


— 1 84 




7.22 


3.96 


1.15 


—2.51 


UUO\J\J .OiJ 


2.01 


0.43 


—0.64 


— i .oo 


UOOOO .Zr 1 


7.86 


4.61 


1.81 


— 1.91 


f\f\^Af\ 49 


2.20 


0.60 


—0.53 


1 7fi 




6.58 


3.26 


0.52 


—3.02 


UU 1 1 o.oo 


1.94 


0.44 


—0.59 


— 1.(0 




6.26 


2.95 


0.20 


—3.31 




1.80 


0.32 


—0.66 


— i .oO 


08734.32 


6.18 


2.85 


0.09 


—3.42 


06828 53 


1.69 


0.24 


—0.74 


— i .oD 


08992.56 


7.77 


4.43 


1.54 


—2.24 


06834 53 


1.65 


0.20 


—0.76 


1 Q9 


09000.47 


7.85 


4.49 


1.60 


—2.17 


06846 50 


1.58 


0.14 


—0.80 


1 Q/1 


09146.20 


8.30 


5.03 


2.13 


—1.68 


06877 45 


1.45 


0.05 


—0.88 


9 

— Z.Uo 


09501.19 


6.60 


3.29 


0.43 


—3.21 


06898 38 


1.39 


-0.02 


—0.92 


—2 10 


09674.61 


7.88 


4.67 


1.78 


—2.11 


06936 27 


1.39 


-0.06 


—0.96 


9 n^^ 

Z.UO 


09820.30 


8.22 


5.07 


2.19 


— 1.68 


06967 28 


1.51 


-0.01 


—0.94 


9 01 
— Z.Ui 


09888.20 


7.49 


4.31 


1.46 


—2.34 


07121 57 


2.21 


0.63 


—0.48 


— 1 .uo 


10111.55 


6.32 


3.11 


0.28 


—3.36 


07149 55 


2.12 


0.53 


—0.47 


1 7^ 
— i . / 


10126.55 


6.34 


3.12 


0.29 


—3.32 


071 76 55 


1.88 


0.38 


—0.62 


1 79 
1 . i z 


10468.50 


8.20 


5.13 


2.25 


— 1.73 


07191 54 


1.83 


0.34 


—0.65 


1 79i. 


10504.41 


8.07 


4.97 


2.09 


— 1.85 


0721 9 49 


1.80 


0.30 


—0.65 


1 7Q 

— -L. 1 O 


10805.60 


6.76 


3.47 


0.57 


—3.18 


07240.41 


1.72 


0.26 


—0.67 


—1.74 












07288 38 


1.73 


0.27 


—0.63 


1 7 A 


IRAS 10131 + 3049 






07507 56 


1.62 


0.14 


—0.81 


— 1 91 


10111.54 


6.73 


3.98 


1.74 


—0.99 




1.78 


0.26 


—0.73 


— 1 .oo 


10468.53 


5.48 


2.85 


U. / D 


1 7S 
— ±. / o 


U / 000. Oo 


1 Q9 


u.oy 


— U.DO 


1 S9 
— 1 .oZ 


10504.43 


5.35 


2.71 


0.65 


-1.84 


07618.38 


2.22 


0.63 


-0.47 


— 1 .uo 


10980.20 


6.90 


4.13 


1.89 


-0.84 


07685.29 


2.19 


0.59 


-0.49 


— 1 62 












07896.56 


1.69 


0.26 


-0.68 


— 1.85 


IRAS 10491 - 2059 






08029.29 


1.81 


0.31 


-0.68 


-1.84 


05008.49 


2.48 


0.85 


-0.33 


-1.59 


08073.21 


2.02 


0.46 


-0.56 


-1.71 


05033.48 


2.29 


0.72 


-0.37 


-1.59 


10483.50 


2.48 


0.88 


-0.20 


-1.62 


06074.56 


2.03 


0.55 


-0.53 


-1.71 


10500.48 


2.42 


0.83 


-0.23 


-1.63 


06111.49 


1.90 


0.48 


-0.58 


-1.75 


10511.43 


2.44 


0.82 


-0.23 


-1.61 


06132.45 


1.81 


0.37 


-0.64 


-1.78 


10520.47 


2.37 


0.77 


-0.28 


-1.65 


06151.43 


1.75 


0.30 


-0.68 


-1.83 


10590.30 


2.20 


0.65 


-0.38 


-1.89 


06181.33 


1.75 


0.33 




-1.81 


10615.23 






-0.42 




06223.30 


1.78 


0.34 




-1.72 


10615.23 


2.27 


0.65 


-0.40 


-1.80 


06270.21 


1.74 


0.31 


-0.67 


-1.76 


10639.21 


2.29 


0.66 


-0.39 


-1.81 


06455.55 


1.59 


0.09 


-0.87 


-2.09 


10639.22 


2.27 


0.62 


-0.45 





DATE 


J 


H 


K 


L 


DATE 


J 


H 


iX 


T 

Li 


IRAS 10491 -2059 .... 


.continued 


IRAS 17049 - 2440 






10802.56 




0.74 


—0.35 




09824 50 




7.66 




1 24 


10802.56 


2.41 


0.73 




—1.67 






8.15 






10809.58 


2.41 


0.73 


— u.oo 


— J..UO 


09911.47 




8.35 


^ QA 
0.O4 


1 7fi 












09945.31 




8.65 


5 fil 


1 8Q 
1 .oy 


IRAS 12447 + 0425 






10180.63 




9.56 


zl^ 
D.40 


z.oo 


nfif^Q7 94 
uuoy t . zi 


4.63 


2.94 


1.71 


0.23 


10208.52 




9.60 


QS 


9 (\0 
Z.DZ 


uuouo. oo 


5.08 


3.42 


2.07 


0.46 


10233.46 




9.46 


fi 95 


2 51 


uuozo.uz 


4.90 


3.28 


1.97 


0.41 


10255.42 




9.18 


O.UO 


9 "IQ 

z.oy 


V\({RAA ^fk 


4.63 


3.04 


1.77 


0.20 


10295.39 




8.63 


5 5*^ 
o.oo 


1 Q1 


U\JO 1 0.10 


4.30 


2.74 


1.52 


—0.01 


10317.29 




8.28 


O.Zo 


1 7(\ 


nfi8Q4 A^ 


4.22 


2.64 


1.45 


—0.10 


10362.23 




7.94 


4 Q5 


1 44 


uuoy y . oy 


4.12 


2.62 


1.43 


—0.08 


10531.60 




7.53 


4 57 


1 1 


uuy'io.ou 


4 15 


2 52 


1.37 


—0.09 


10692.28 




8.32 


^ QA 
0.O4 


1 

i .OU 


uuyu t . oo 


4.22 


2.58 


1.41 


—0.08 


10912.57 


12.56 


9.32 


u.ou 


9 7*^ 
Z. ( o 


071 79 fiO 
U / 1 / z.uu 


5 97 


'\ 54 


2.22 


0.66 


10980.47 




9.29 


R Q7 


9 71 
Z. ( 1 


OVI 01 fiO 

u 1 iy i.DU 


^ nn 
o.uu 


o.o4 


2.07 


0.47 












07236.45 


4.95 


3.26 


2.01 


0.43 


IRAS 17119 + 0859 






07263.41 


4.90 


3.22 


2.00 


0.43 


09912.42 


4.51 


3.16 


2 25 


n 81 

U.O -L 


07368.19 


3.99 


2.35 


1.27 


—0.22 


10214.54 


5.38 


3.89 


2.65 


0.95 


07539.56 


5.60 


3.74 


2.29 


0.71 


10237.45 


4.79 


3.44 


2 27 


n R9 

u.uz 


07605.56 


5.56 


3.75 


9 Qfl 

z.ou 


n ^^^^ 












07629.43 


5.45 


3.63 


2.20 


0.59 


IRAS 17297 + 1747 






07669.27 


5.28 


3.48 


2.10 


0.56 


09912.41 


9.87 


6.84 


A Pin 

4.0U 


i .^0 


07690.29 


5.31 


3.51 


2.11 


0.53 


10214.55 


8.30 


5.27 


9 09 

z.yz 


U.Uo 


07903.59 


4.78 


2.97 


1.69 


0.17 


10237.47 


8.26 


5.21 


2.88 


0.08 


07989.43 


5.62 


3.73 


2.29 


0.73 












08299.56 


4.41 


2.61 


1.40 


-0.15 


IRAS 17360 - 3012 






08323.53 


4.63 


2.78 


1.53 


-0.06 


09824.52 


9.58 


6.36 


4.36 


2.12 


08388.36 


5.31 


3.35 


1.99 


0.41 


10210.57 


10.50 


7.28 


5.03 


2.36 


08457.21 


5.76 


3.81 


2.35 


0.74 


10238.51 


10.24 


7.05 


4.84 


2.20 


08735.45 


4.51 


2.66 


1.36 


-0.20 


10358.25 


8.73 


5.78 


3.76 


1.37 


09172.26 


4.70 


2.69 


1.34 


-0.20 


10590.56 


7.87 


5.06 


3.16 


0.93 


09888.22 


4.99 


3.22 


1.97 


0.44 


10618.48 


7.83 


5.02 


3.14 


0.95 


10112.50 


4.59 


2.74 


1.51 


0.09 


10719.26 


7.97 


5.12 


3.27 


1.07 


10504.52 


4.15 


2.40 


1.27 


-0.19 


10976.49 


9.38 


6.41 


4.48 


2.13 



b 

e 

03 
I 

03 

3 

01 

>^ 
ta 

03 



Co 
o 



5- 
O 

o 
o 



Table 2. ...continued. Near-infrared data. The Date is given as JD-2440000.00. 
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CO 
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o 
o 

i 
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From the number of references listed in Tables 1, CW 
Leo stands out as the best studied source in this sample. 
About 50 molecular species have been detected in its en- 
velope (Wallerstein & Knapp 1998) which is more than for 
any other source. The full sample consists mainly of cool 
variable stars with a few possible post-AGB objects. 



3 OBSERVATIONS 

3.1 Near-infrared Photometry 

Near-infrared observations were obtained in the JHKL 
bands (1.2, 1.65, 2.2 and 3.45^im) for 42 of the 45 south- 
ern sources, defined henceforth as stars with S < -|-30°. 
Some of these 42 stars were selected exclusively from JK89, 
while the others were already part of other monitoring pro- 
grammes. The remaining three southern stars (06500-1-0829, 
17513 - 2313 and RAFGL 1406) were not observed because 
they were in crowded fields and confusion rendered observa- 
tions with a single channel photometer impossible (note that 
Winfrey ct al. (1994) suggest that 17513 - 2313 is actually 
a supergiant and it is therefore not referred to again). 

The near-infrared observations for the 42 southern 
sources, accurate to ±0.1 mag and better (0.03 mag in JHK 
and 0.05 mag in L in the majority of cases), are listed in Ta- 
ble 2. The magnitudes from the 0.75 m reflector arc on the 
SAAO system as defined by Carter (1990). The HKL bands 
for the 1.9 m reflector are assumed to be identical to those 
of the 0.75 m reflector, while the J-magnitude is converted 
to the SAAO system. 

In some instances, no observations arc listed because the 
sources were too faint to be reliably mcEisured through the 
relevant filter with a particular telescope. The times of the 
observations are given as Julian dates from which 244 0000 
days has been subtracted. In addition to the data from Ta- 
ble 2, published SAAO observations from Whitelock et al. 
(1994, 1997a) and Lloyd Evans (1997) have been used in 
the analysis, together with several unpublished observations 
from Lloyd Evans (private communication) with a sum total 
of 848 measurements. 

Mean magnitudes calculated from the average of ob- 
served values at maximum and minimum light are listed in 
Table 3 and compared with the Fourier mean in appendix A. 
These two averages compare fairly well, differences were less 
than 0.3 mag for the J-band and less than about 0.1 mag for 
the HKL bands. Table 3 also lists mean is'-magnitudes cal- 
culated as above from values listed in the literature, for the 
20 remaining stars with no near-infraxed SAAO photometry. 



3.2.2 Low Resolution Spectra (LRS) 

The IRAS Low Resolution spectral classifications, from Ol- 
non &: Raimond (1986) or Loup et al. (1993), are listed in 
Table 4. Note that 18135 — 1641 does not have a published 
classification so we examined the spectrum in the IRAS Low 
Resolution Spectra electronic database at the University of 
Calgary. It shows a weak 10/xm emission feature, but no 
IS/um feature, which is present in the spectra of all the other 
oxygen-rich stars. It is actually a supergiant as discussed be- 
low (section 4). 



4 VARIABILITY 

Using the data in Table 2, the southern sources were classi- 
fied as variable in the near-infrared if the difference between 
the maxinmm and mininmm magnitudes at K and/or at 
L was greater than 0.1 mag; thirty seven of the stars axe 
variable and 2 non- variable. The only two decisively non- 
variable stars, 18135 - 1641 and 18204 - 1344, are the super- 
giants (see section 2). There were insufficient observations 
to determine if the other 6 stars arc variable or not. Table 

5 lists the variability information. 

Period determinations, using Fourier transforms of the 
K light curve, were made for all variables with measure- 
ments on at least 10 epochs. Plausible, although not defini- 
tive, periods were also established for some sources with as 
few as 8 measurements. For the very red source 17411 — 3154, 
a period could only be estimated from the L data. Periods, 
listed in Table 5, were thus determined for a total of 20 stars, 
and phased K light curves for these are shown in Fig. 1. 
Where possible the Fourier transform analysis was carried 
out at the other wavelengths and in all cases the periods 
found were within 2 percent of those listed in Table 5. No 
period determination was possible for 19126 — 0708 (W Aql), 
although it had been observed on 12 epochs. This is because 

6 of the epochs were between 1968 and 1982 when this car- 
bon star was undergoing a dust obscuration episode (Mattei 
& Foster 1997). Such an event would disrupt the apparent 
periodicity. 

The average of periods found in the literature are listed 
in Table 5, with references, for the 47 stars concerned. Where 
there are periods in common none of them differs signifi- 
cantly from those determined here. In the following analysis 
our own determinations, where available, are used in pref- 
erence to others, except for 17360 - 3012 and 18348 - 0526 
where our determination is uncertain due to the small num- 
ber of measurements available. 



3.2 Mid- and Fcir-infrared Observations 

3.2.1 Photometry 

Infrared fluxes at 12, 25 and 60^m were taken from the IRAS 
Point Source Catalog (PSC) and colour corrected as speci- 
fied in the IRAS Explanatory Supplement (ES). The results 
are listed in Table 4 together with the IRAS fiux-quality 
designations and the IRAS variability index (V.L), which 
gives the probability that the source is variable in units of 
0.1 (see IRAS ES). Table 4 also contains ll/nm magnitudes 
from the Revised Airforce Geophysics Laboratory, RAFGL, 
survey. 



4.1 Phased Light-curves 

Using the periods found for these 20 stars, the relative 
data were phased, arbitrarily assuming zero phase at JD 
2440000.00, where JD is the heliocentric Julian date. These 
phased light-curves are plotted in Fig. 1 together with least 
squares fits to the data of the form: 
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Table 3. Average neax-infraxed magnitudes 
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Table 4. Mid-infrared and IRAS Data. 
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RAFGL 1406 -3.0 
RAFGL 2688 -2.6 

All stars have LRS classifications taJjen from IRAS Atlas of Low Resolution Spectra 

(Olnon & Raimond, 1986), except in the cases marked which were taJsen from Loup et al. (1993) 

f-Flux Density Uncertainty, as explained in the IRAS Explanatory Supplement. 
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Table 5. Variability data 



NAME 



P (days) 



Piit (days) VAR AJ AH AK AL 



IRAS 00042 + 
IRAS 01037 + 
IRAS 01159 + 
IRAS 02270 - 
IRAS 02316 + 
IRAS 02351 - 
IRAS 03229 + 
IRAS 03507 + 
IRAS 04307 + 
IRAS 04566 + 
IRAS 05073 + 
IRAS 05411 + 
IRAS 05559 + 
IRAS 06176 - 
IRAS 06300 + 
IRAS 06500 + 
IRAS 08088 - 
IRAS 09116 - 
IRAS 09429 - 
IRAS 09452 + 
IRAS 10131 + 
IRAS 10491 - 
IRAS 12447 + 
IRAS 17049 - 
IRAS 17119 + 
IRAS 17297 + 
IRAS 17360 - 
IRAS 17411 - 
IRAS 18009 - 
IRAS 18040 - 
IRAS 18194 - 
IRAS 18240 + 
IRAS 18333 + 
IRAS 18348 - 
IRAS 18349 + 
IRAS 18397 + 
IRAS 18398 - 
IRAS 18413 + 
IRAS 18560 - 
IRAS 19008 + 
IRAS 19059 - 
IRAS 19093 - 
IRAS 19126 - 
IRAS 19175 - 
IRAS 19321 + 
IRAS 20077 - 
IRAS 20396 + 
IRAS 20440 - 
IRAS 20570 + 
IRAS 21032 - 
IRAS 21286 + 
IRAS 21320 + 
IRAS 21456 + 
IRAS 23166 + 
IRAS 23320 + 
IRAS 23496 + 
RAFGL 1406 
RAFGL 2688 



1219 
1219 
7220 
2619 
6455 
2711 
4721 
1115 
6210 
5606 
5248 
6957 
7430 
1036 
6058 
0829 
3243 
2439 
2148 
1330 
3049 
2059 
0425 
2440 
0859 
1747 
3012 
3154 
2019 
0941 
2708 
2326 
0533 
0526 
1023 
1738 
0220 
1354 
2954 
0726 
2219 
3256 
0708 
0807 
2757 
0625 
4757 
0105 
2714 
0024 
1055 
3850 
6422 
1655 
4316 
6131 



645 
389 
480 
469 



570 
670 
640 
652 

530 
441 
775 



(1120) 
1440 



690 

795 
(1500) 



600 

575 

380 
675 



750-^ 

6521.3,5,6 

612^ 
3872.4,5 

534'^ 
4751,3 

5375.6 
4541,3,5,6 

5575 
629^.6 

5235 



5275 
5712 

6433,6 

6352.S.' 

640^ 

5302.5 
4335 



520^ 
1150^ 



15703.5,6 
500^ 
511^ 

590^ 

5772 

511^.5 

3725 

4905 

6762 

6256 

660^.5 

42l5 

750" 

4545 

4575 

4705 

6225 

7OOI.2 

620^ 

4196 



Y 
Y 
Y 
Y 



Y 
Y 
Y 
Y 
Y 
Y 
Y 
Y 
Y 
Y 
Y 
Y 

Y 
Y 
Y 
Y 
Y 
Y 
Y 
Y 

Y 
Y 
Y 
Y 
Y 
Y 
Y 
Y 

Y 
Y 
Y 
Y 



2.87 2.15 1.69 1.34 

1.06 0.88 0.68 0.55 

1.56 1.24 0.92 0.70 

1.96 1.45 1.08 0.94 



2.09 



2.50 



1.74 1.47 1.19 



2.02 2.03 1.54 
1.82 1.56 1.39 1.35 



2.14 2.20 2.11 



1.81 



1.40 1.22 0.99 0.91 
2.25 2.13 1.81 



1.90 1.99 1.77 1.55 
1.11 



1.66 1.57 1.35 1.20 

3.00 2.08 1.48 1.17 
3.19 2.11 



1.25 1.09 .93 0.75 
2.32 1.81 1.35 1.00 

1.21 0.98 .78 0.73 



1.45 1.22 



REFERENCES-(l) Whitclock ct al. (1994): (2) Le Bertre (1993a); 
P.N. (1985)(GCVS); (6) Jones ct al. (1990). 



(3) Le Bertre(1993b) ; (4) Whitelock et al. (1997a); (5) Kholopov, 



© 0000 RAS, MNRAS 000, 000-000 



12 E. A. Olivier et al. 



Figure 1. Phased Light-curves 
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Figure 1. Phased Lightcurves continued. 
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M = ^-sin((/> - 4,0) + Mav (1) 

where M is the magnitude in a particular band, AM the 
peak to peak amphtude, 4> related to the phase, p, hy (j) = 
2-np and cjio the phase angle at zero phase. Mav is the mean 
magnitude. 

One immediate observation from most of the phased 
light-curves is that at a given phase there is a range of mag- 
nitudes, rather than one sharply defined magnitude. This is 
not, in general, due to observational error, but is contributed 
to by various effects. First, Miras and OH/IR stars are er- 
ratic pulsators in the sense that their light curves do not ex- 
actly repeat from one cycle to the next at visual or infrared 
wavelengths (Koen & Lombard 1995; Whitelock, Marang & 
Feast 2000) Secondly, some stars, particularly carbon stars 
go through faint phases from time to time due to variable 
dust obscuration, e.g. 02270-2619 (R For, see Whitelock et 
al 1997a). Thirdly, a very small number of stars show tem- 
poral or long term changes in luminosity, amplitude and/or 
period, possibly due to evolutionary effects (Whitelock 1999 
and references therein). 

Many of the light curves in Fig. 1 show clear departures 
from simple sinusoids and would be better fitted by the in- 
clusion of one or two harmonics, in particular, 09429 — 2148 
and 12447 + 0425. The most strikingly odd curve is that of 
10491 - 2059 (V Hya) which is discussed below. 

4.1.1 V Hya 

V Hya is a well studied carbon star and semi-regular vari- 
able. It has a bipolar circumstellar outfiow, and models in- 
voke rotation of the central star or the presence of a binary 
companion to explain bipolar outflows (Barnbaum, Morris 
& Kahane 1995). 

Fig. 2a illustrates the light curve for V Hya using K 
magnitudes measured from SAAO, including material from 
Lloyd Evans (1997 and private communication) together 
with the data from Table 2. In addition to the regular pulsa- 
tion V Hya experienced two faint phases around JD 244 3500 
and 244 9500. A Fourier analysis of the K data when the 
star is bright, between JD 244 6000 and 244 8250 provide a 
period of 530 days, identical to the value given in the GOVS. 
Fig. 2b shows the data phased at this period and illustrates 
the double peak caused by harmonics in the light curve. 

Knapp et al. (1999) analysed 3260 AAVSO visual mag- 
nitude estimates covering almost 35 years. They confirm 
and refine the very long period found by Mayall (1965) at 
6160 ± 400 days. The two faint phases seen in Fig. 2a are 
consistent with this long period. The amplitudes of these 
6160 day variations are at least 2.4, 2.1, 1.7 and 0.7 mag, 
at JHKL, respectively, while AV ~ 3.5 mag (Knapp et al. 
1999). 

Barnbaum et al. (1995) found evidence for rapid rota- 
tion in V Hya, using spectral line broadening analysis. They 
suggest that the long term variation can be attributed to 
a coupling between the radial pulsation, with period of 530 
days, and the rotation of the star with a very similar period. 
It is not clear how this interaction would give rise to the very 
long period without strongly modulating the shorter period, 
or why it should result in the observed reddening of the star 
during the faint phases. 

Lloyd Evans (1997) suggested that the long period may 
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Figure 2. (a) K light curve of V Hya. (b) data between 
JD 244 6000 and 244 8250 phased at 530 days. 

be the result of obscuration events like those observed in R 
For. Although the colours would be consistent with this, the 
regularity of the long period (Knapp et al. 1999) is quite un- 
like the effect observed in R For and this explanation there- 
fore seems unlikely. 

Knapp et al. (1999) on the other hand, suggest that 
the long-period dimming of this star is due to a thick dust 
cloud orbiting the star and attached to a binary compan- 
ion. Indeed, there is evidence for large numbers of late-type 
binaries in the LMC (Wood et al. 1999), so this may not 
be rare occurrence. The colour variations of V Hya are less 
extreme than those expected from optically thin interstellar 
dust, but they could be explained by large particles. 

4.2 Period Distributions 

The period distributions for the oxygen- and carbon-rich 
stars are plotted in Fig. 3. Most of the stars in the sam- 
ple have periods between 500 and 700 days and almost all 
have periods between 300 and 900 days. There are only three 
extreme OH/IR stars with periods in excess of 1000 days. 
Fig. 3 can be compared with Whitelock et al. (1994) fig. 8 
which shows the period distribution of IRAS selected Miras 
in the South Galactic Cap. Most of the Cap sample have 
150 < P < 450 days, i.e. shorter than the sample under 
discussion, but longer than optically selected Miras. Longer 
period Miras have more massive progenitors and are there- 
fore lower scale heights. It is not surprising to find more of 
them in the solar neighbourhood than in the Galactic Cap. 

A comparison of the period distributions of the oxygen- 
and carbon-rich stars is instructive. There are no carbon- 
rich stars with periods longer than 1000 days. From Fig. 3 
it can be seen that the maximum of the carbon-star period 
distribution seems to be slightly shifted to longer periods 
than that of the oxygen-rich stars. A Kolmogorov-Smirnov 
(KS) test, when applied to these two distributions for pe- 
riods less than 1000 days, yields a probability that the two 
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Figure 3. Period distributions 

distributions come from the same parent population of 0.40. 
This probability (> 0.05) implies that there is no significant 
difference in the period distributions for the carbon- and 
oxygen- rich stars. 

The KS test cannot be used if the very long period stars 
are included since it is not very sensitive to differences in the 
tails of the distributions, as is well known. In this case the 
F-test and then the appropriate version of the student t-test 
is applied. The F-test gives a probability that the variances 
of the two distributions are the same of 2.5 x 10~*, i.e. very 
unlikely, which is obvious from Fig. 3. The appropriate ver- 
sion of the t-test, which takes differences in variances into 
account, gives a probability that the two distributions have 
the same mean of 0.24. This probability (> 0.05) implies 
that there is no statistical evidence for any difference in the 
mean of the two distributions. 

4.3 Amplitudes 

Table 5 lists the stars for which periods were determined 
together with their amplitudes, from equation 1, at JHKL. 
These amplitudes are plotted against period in Fig. 4, where 
a clear correlation can be seen at all wavelengths. 

The amplitude clearly decreases with wavelength (see 
also Feast et al. 1982). Quantitatively, the average ampli- 
tudes, relative to that at J are, for H, K and L respec- 
tively, 0.84 ± 0.1, 0.69 ± 0.14 and 0.59 ± 0.14. The linear 
correlation coefficient between the average relative ampli- 
tude and wavelength is —0.93 with a correlation probability 
of Prob = 0.07 (using the relative amplitude at J of 1). This 
probability is considered as marginally significant. The KS 
test gives the probability for the null hypothesis that the 
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Figure 4. Light-Amplitude vs. pulsation period, the lines show 
least squares fits; filled and open circles represent carbon- and 
oxygen-rich stars, respectively. 

light amplitude distributions for the carbon- and oxygen- 
rich stars (with P < 1000 days) came from the same parent 
population at J, H, K and L of 0.25, 0.92, 0.51 and 0.52 
respectively. Hence there is no statistical reason to suspect 
that these carbon- and oxygen-rich stars have significantly 
different light-amplitudes in the near-infrared. 



5 COLOURS 

5.1 Average Colours 

Average near-infraxed colours, i.e. consisting of JHKL mag- 
nitudes only, were calculated for each star using the magni- 
tudes at the dates of minimum and maximum /f- light. This 
average is compared to a Fourier-fit colour in Appendix A. 
The two values are quite comparable with mean differences 
for all the near-infrared colours not exceeding 0.1 mag. Note 
that due to the very small phase coverage of 17411 — 3154 
in the ii'-band, the derived K — L colour of this star may 
not be a good representative of the average K — L colour. 

The IRAS colours were obtained using the data in Ta- 
ble 4, where the colour-corrected IRAS fluxes were converted 
to magnitudes using the zero points listed in the IRAS-ES. 
Colours consisting of JHKL and IRAS photometry, were 
calculated by combining average magnitudes in the JHKL 
bands (Table 3) with the IRAS magnitudes. 

For the northern sources and the three non-IRAS 
sources the average K magnitudes were taken from Table 3, 
to calculate the K — [12] colour. For the non-IRAS (caxbon- 
rich) sources the average 12/xm magnitudes were estimated 
from the average [11] — [12] = 0.1±0.3 colour for the carbon- 
rich sources and the 11 fim magnitudes listed in Table 4. All 
the derived infrared colours are given in Table 6. 
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Table 6. Infrared Colours 



NAME 


J - H 


H — K 


J - K 


K — L 


K - [12] 


[12] - [25] 


[25] - [60] 


IRAS 00042 + 4248 










5.37 


1.06 


0.04 


IRAS 01037 + 1219 


3.51 


2.26 


5.77 


2.48 


6.26 


1.21 


0.25 


IRAS 01159 + 7220 










4.65 


0.96 


-0.22 


IRAS 02270 — 2619 


1.71 


1.33 


3.04 


1.47 


3.54 


0.19 


0.24 


IRAS 02316 + 6455 










5.87 


1.07 


-0.16 


IRAS 02351 - 2711 


1.25 


0.69 


1.94 


0.79 


3.66 


0.86 


-0.32 


IRAS 03229 + 4721 










4.60 


0.42 


0.18 


IRAS 03507 + 1115 


1.96 


1.16 


3.12 


1.32 


4.81 


0.68 


—0.27 


IRAS 04307 + 6210 










3.90 


0.42 


0.12 


IRAS 04566 + 5606 










3.69 


0.45 


0.23 


IRAS 05073 + 5248 










5.22 


1.56 


0.41 


IRAS 05411 + 6957 










4.30 


0.84 


-0.27 


IRAS 05559 + 7430 










3.36 


0.79 


-0.03 


IRAS 06176 - 1036 


1.66 


1.58 


3.24 


2.11 


6.43 


1.40 


0.78 


IRAS 06300 + 6058 










3.48 


1.06 


0.22 


IRAS 06500 + 0829 










4.34 


0.82 


0.54 


IRAS 08088 — 3243 


2.59 


2.12 


4.71 


2.62 


6.37 


0.60 


0.18 


IRAS 09116 - 2439 




3.75 




3.56 


9.40 


0.78 


0.20 


IRAS 09429 — 2148 


1.86 


1.34 


3.19 


1.62 


5.39 


1.30 


—0.16 


IRAS 09452 + 1330 


3.20 


2.82 


6.02 


3.74 


9.09 


0.64 


0.37 


IRAS 10131 + 3049 


2.70 


2.15 


4.85 


2.61 


6.24 


0.38 


0.29 


IRAS 10491 - 2059 


1.79 


1.27 


3.06 


1.55 


3.69 


0.51 


0.24 


IRAS 12447 + 0425 


1.85 


1.29 


3.15 


1.54 


3.80 


0.22 


0.17 


IRAS 17049 - 2440 




3.04 




3.62 


9.02 


0.90 


0.33 


IRAS 17119 + 0859 


1.42 


1.08 


2.50 


1.57 


5.10 


1.18 


-0.28 


IRAS 17297 + 1747 


3.04 


2.33 


5.37 


2.93 


6.73 


1.13 


0.04 


IRAS 17360 - 3012 


3.01 


2.07 


5.08 


2.43 


6.28 


1.79 


0.26 


IRAS 17411 - 3154 








6.25 


13.99 


2.14 


1.04 


IRAS 18009 — 2019 


1.73 


1.02 


2.74 


1.01 


3.78 


1.11 


-0.07 


IRAS 18040 — 0941 


2.46 


1.74 


4.20 


1.95 


4.43 


0.45 


0.36 


IRAS 18194 — 2708 


3.20 


2.39 


5.59 


2.80 


7.06 


0.40 


0.45 


IRAS 18240 + 2326 




3.24 




3.93 


9.21 


0.92 


0.14 


IRAS 18333 + 0533 


4.26 


2.48 


6.74 


2.44 


6.25 


1.44 


0.32 


IRAS 18348 - 0526 








6.27 


11.19 


1.94 


1.44 


IRAS 18349 + 1023 


1.55 


0.90 


2.45 


1.09 


3.98 


0.58 


0.20 


IRAS 18397 + 1738 


2.14 


1.74 


3.88 


2.12 


4.92 


0.56 


0.40 


IRAS 18398 - 0220 


2.40 


1.72 


4.12 


1.98 


4.83 


0.57 


0.23 


IRAS 18413 + 1354 


1.81 


0.97 


2.78 


0.98 


3.97 


1.13 


—0.19 


IRAS 18560 - 2954 


1.47 


0.84 


2.31 


1.07 


3.99 


0.76 


0.12 


IRAS 19008 + 0726 


2.63 


2.04 


4.67 


2.45 


5.38 


0.49 


0.20 


IRAS 19059 - 2219 


1.72 


1.04 


2.77 


1.12 


4.45 


1.15 


-0.07 


IRAS 19093 — 3256 


1.31 


0.70 


2.01 


0.76 


3.61 


1.02 


-0.08 


IRAS 19126 - 0708 


1.31 


0.79 


2.10 




4.61 


0.60 


-0.02 


IRAS 19175 - 0807 


2.70 


2.04 


4.74 


2.46 


5.46 


0.67 


0.38 


IRAS 19321 + 2757 


2.42 


1.89 


4.31 


2.29 


5.06 


0.73 


0.30 


IRAS 20077 - 0625 


2.89 


1.90 


4.79 


1.98 


6.13 


1.24 


0.16 


IRAS 20396 + 4757 










3.60 


0.35 


0.22 


IRAS 20440 - 0105 


1.18 


0.56 


1.74 


0.65 


3.16 


0.89 


-0.14 


IRAS 20570 + 2714 


2.96 


2.43 


5.39 


2.88 


5.83 


0.76 


0.22 


IRAS 21032 - 0024 


1.92 


1.35 


3.27 


1.51 


3.71 


0.44 


0.14 


IRAS 21286 + 1055 


1.38 


0.86 


2.24 


1.01 


3.48 


1.04 


0.00 


IRAS 21320 + 3850 










3.75 


0.38 


0.26 


IRAS 21456 + 6422 










3.17 


0.99 


-0.10 


IRAS 23166 + 1655 








6.23 


14.05 


1.43 


0.61 


IRAS 23320 + 4316 










7.23 


0.66 


0.34 


IRAS 23496 + 6131 










4.98 


1.08 


0.03 



RAFGL 1406 5.1 
RAFGL 2688 11.1 
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Figure 5. J — H versus H — K; filled circles represent carbon- 
rich stars, open circles oxygen-rich stars and five-point star S- 
typc stars. The solid line represent loci of blackbodies of different 
temperatures. 
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Figure 6. J-K versus K-L. Symbols as in Fig. 5. 



5.2 Colour-Colour Diagrams 

Infrared two colour diagrams are illustrated in Figs. 5 to 8. 
There is a clear separation of the oxygen- and carbon-rich 
stars in these diagrams with the carbon stars falling to the 
lower right of the oxygen-rich stars in the near-infrared two 
colour diagrams (Figs. 5 and 6) and to the left of them in 
the combined diagram, Fig. 7. The JHKL colours of these 
AGB stars are more extreme than those of the well studied 
Miras (e.g. Feast et al. 1982), or even of IRAS selected Miras 
in the Galactic Cap (Whitelock et al. 1994). 



5.2.1 The K - [12] versus [12] - [25] Diagram 

Of particular interest is the K — [12] vs. [12] — [25] diagram 
(Fig. 7). The colours of the oxygen- and carbon-rich stars 
separated rather clearly in this diagram. This has been noted 
previously (Le Bertre et al. 1994) and shown to bo the result 
of the difference in the ratios of the near-infrared to mid- 
infrared opacities for the carbon- and oxygen-rich dust. 

The single carbon star which falls among the oxygen- 
rich stars is 06176—1026 (the Red Rectangle), a well studied, 
but peculiar star. It has an IRAS LRS classification of 80 
indicating the presence of the 11.3/im line in emission in 
contrast to the other carbon stars which show 11.3/im ab- 
sorption, hence its position in Fig. 7. It is thought to be a 




[12]-[25] 

Figure 7. K-[12] versus [12]-[25]. Symbols as in Fig. 5. Dotted 
line separates the chemical types and is given by equation 3. 

pre-planetary nebula with an oxygen-rich circumstellar disk 
(see section 2.1). 

The dashed line in Fig. 7 has been positioned so as to 
separate the carbon- and oxygen-rich stars as far as possible, 
it is described by: 

{K - [12]) = 11.5([12] - [25]) - 3.8 (2) 
Thus if we define 5c /o as follows: 

6c/o = 11.5([12] - [25]) -{K- [12]) - 3.8, (3) 

The two chemical types can be separated by a single 
parameter, 5c/0: such that for carbon stars 5c/o < ^ and 
for most oxygen- rich stars 5c /o > 0. This parameter might 
be useful for separating these objects over the colour range 
discussed. On this sample, obviously the one used to define 
the parameter, only 6 percent of stars were misclassified. 

The two S-type stars were omitted from the above dis- 
cussion, we note that one of them falls with the carbon- and 
the other with the oxygen-rich stars. 

5.2.2 IRAS Two Colour Diagram 

Van der Veen & Habing (1988) demonstrated that the posi- 
tion of an AGB star in an IRAS two-colour diagram, of the 
type shown in Fig. 8, was indicative of its chemical type and 
the thickness of its dust shell. Combined with information 
on variability and IRAS spectral-type the diagram provides 
a useful diagnostic. The colour indices for this plot (C12/25, 
C26/60) are defined by: 

C25/60 = 2.51og(f5il!n), 

where the IRAS fluxes are not colour corrected. 

Following van der Veen & Habing the diagram is divided 
into regions, each of which is dominated by a certain type of 
object. Stars with low mass-loss rates will fall in the lower 
left of the diagram and increasingly higher mass-loss rates 
are found as one moves to the upper right. 

Most of the oxygen-rich stars under discussion occupy 
region Ilia which van der Veen & Habing found to con- 
tain variable stars with moderate oxygen-rich dust-shells. 



© 0000 RAS, MNRAS 000, 000-000 



18 E. A. Olivier et al. 




Figure 8. C25/6O versus Ci2/25- Symbols as in Fig. 5. 18135 — 
1641 has only an upper limit to its C25/60 colour as indicated by 
the arrow. 



Four stars are located in regions Illb and IV, indicating 
thick to very thick oxygen-rich dust-shells and high mass- 
loss rates. Four oxygen-rich stars axe located in region VII, 
an area dominated by carbon stars, but where van der Veen 
& Habing also found oxygon-rich stars. 

Most of the LRS spectra for oxygen-rich stars in this 
sample are classified as 2n, indicating that the lO/itm silicate 
feature is in emission. The two oxygen-rich variables with the 
reddest [12] - [25] colour, 17411 - 3154 and 18348 - 0526, 
have LRS classifications of 3n, i.e. their spectra show the 
10/im silicate feature in absorption. In these stars the dust- 
shell has become so cool and optically thick that the silicate 
feature has gone into absorption. 

Most of the carbon-stars fall in region VII, which is 
characterised by variable stars with well developed carbon- 
rich dust-shells. All of these stars have an LRS classification 
of 4n, indicating that the SiC 11/im feature is in emission. 
The Red Rectangle, 06176 — 1036, (discussed above) is lo- 
cated in region Ilia, as is 23166 + 1655. The latter object has 
an LRS classification of 02, probably due to self absorption 
of the SiC feature in the very think shell (Loup et al. 1993). 



6 BOLOMETRIC MAGNITUDES AND 
DISTANCES 

6.1 Apparent Bolometric Magnitudes 

Given that these sources were selected on the basis of their 
high fluxes and therefore small distances, together with the 
fact that the bulk of their energy is emitted at infrared 
wavelengths it is reasonable to assume that the effects of 
interstellar extinction are negligible. The apparent bolomet- 
ric magnitude, mtoi, were calculated by integrating under a 
spline curve (Hill 1982), fitted to the J, H, K, L, 12- and 25- 
Hm fiuxes as a function of frequency. The end-points were 
dealt with by extrapolating a line, which joins the J flux 
and the point lying halfway between the H and K fluxes, 
to zero flux. Zero flux at zero frequency was assumed. The 
derived values of ruboi are listed in Table 7 for the sources 
with sufficient photometry. 

For the very red sources 23166-1-1655 and 18348-0526, 
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Figure 9. BC12 versus K — [12]. Symbols as in Fig. 5. Curves 
are those given in equations 4 and 4. 



which have no measurements in the J and/or H bands, zero 
flux at J and H was assumed when making the spline fit. 

To estimate an upper-limit to the amount of fiux not 
accounted for by ignoring the tail in the flux distribution, 
a Wien approximation, r-jAi'^e~'"' (where A and b are con- 
stants), was assumed for flux at frequencies higher than J 
with the same slope as the spline at J. In more than 75 per- 
cent of the cases the difference between the fluxes calculated 
in this way and the simple spline fit was less than 1 percent. 
20440-0105, which has the smallest K-[12] = 3.16, shows 
the maximum difference of 5 percent. Neglecting H2O ab- 
sorption between the J,H,K and L bands (Robertson & 
Feast 1981) will have a negligible effect on the bolometric 
magnitude estimated for stars whose flux originates largely 
from the circumstellar dust-shell. 

The bolometric corrections at 12 /im, BCrz, is plotted 
against the K — [12] colour in Fig. 9, together with curves 
given by: 

BC7i2 = -15.4e-°-^°(^-li^»+8.10, 



and 
BCv 



-7.59e-°«(^-li^»+8.41, 



(4) 



(5) 



which are least-square fits of this form to the oxygen- and 
carbon-rich star data respectively. The rms difference be- 
tween the data and these curves is 0.04 and 0.06 mag re- 
spectively. For the remaining sources ruboi was calculated 
from [12] and K — [12] using equations 4 and 5. 

The apparent luminosities derived here arc typically 
rather smaller than those estimated by JK89, because of 
the colour corrections applied to the IRAS photometry. 



6.2 Absolute Bolometric Magnitudes and 
Luminosities 

For the 47 stars with periods the absolute bolometric magni- 
tudes, Mboi, was calculated from the period- luminosity (PL) 
relation, equation 6, derived for LMC variables (Feast et 
al. 1989) assuming an LMC distance modulus of 18.60 mag 
(Whitelock et al. 1997b) and the results given in Table 7. 



-2.34 log P + 1.26 



(6) 



The average luminosity, for all stars with periods less 
than 1000 days is (0.98±0.17) x 10* Lq, or Mm = -5.23± 
0.19. For stars with unknown periods this mean value is 
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assumed; which seems reasonable since the carbon stars 
seem to have periods less than 1000 days and the remain- 
ing oxygen-rich stars have [12] — [25] colours suggesting that 
they too have periods less than a 1000 days. This value of 
the mean luminosity suggests that the 10* Lq often assumed 
for the luminosity of the tip of the AGB is consistent with 
the application of the PL relation to this type of AGB star. 
Note this value may not be applicable to the Red Rectan- 
gle and the Egg Nebula, since they may have left the AGB 
already. 

It is not really clear at this stage how well the various 
PL relations fit long period variables or if different relations 
should be used for the oxygen- and carbon-rich stars. Pro- 
vided this caveat is born in mind when the results are inter- 
preted it seems reasonable to use the luminosities derived as 
described above. 

Note, in particular, that the PL relation used above 
was derived for stars with periods less than 420 days and 
it may not be appropriate to use it for the present sample. 
There certainly are long period AGB {P > 420 days) stars 
in the LMC that lie above the PL relation (e.g. Feast et 
al. 1989). However, recent work on dust-enshrouded stars in 
the LMC suggest that these AGB stars are over-luminous 
because they are undergoing hot bottom burning and that 
more normal long period {P > 400 day) carbon- and oxygen- 
rich AGB stars fall very close to the PL derived for shorter 
period oxygen- rich stars (Whitelock & Feast 2000). 

6.3 Distances 

Distances to all the stars in the sample (excluding the Red 
Rectangle and the Egg Nebula) were calculated from rriboi 
and Mboi and are listed in Table 7. These distances are typ- 
ically larger than those estimated by JK89. As discussed 
above the apparent luminosities used here are smaller than 
JK89's, while the absolute luminosities are on average the 
same, except for the 3 stars with P > 1000 where they are 
brighter. Thus we find sources with distances up to 2 kpc. 

The other plausible way of estimating distances to 
OH/IR stars is the phase lag method, although it is no- 
toriously difficult and time consuming. Van Langevelde, van 
der Heiden & Schooneveld (1990) measured phcise lag dis- 
tances for two of the stars under discussion. For 01037-1-1219 
(WX Psc) they get 0.74 ± 0.15 kpc, identical to the PL dis- 
tance given in Table 7. For 18348-0526 (OH26.5-f 0.6) they 
give two values, 1.44 ±0.27 and 1.30 ±0.35 kpc, both some- 
what closer than the 1.9 kpc listed in Table 7. Given the 
uncertainties in phase lag and the PL method this differ- 
ence is not problematic. West (1998) has measured a phase 
lag distance for 17411-3154 (OH357.3-1.3) of 1.2±0.4 kpc, 
in agreement with the distance of 0.99 kpc listed in Table 7. 



7 MASS-LOSS RATES 

Mass-loss rates, M, were estimated using the expression 
given by Jura (1987): 



M = 1.7 X 10 "^visrlpcLi '''V,,6oAi„, 



(7) 



where wis is the outflow velocity of the gas in units of 15 
kms^^, rkpc the distance to the star in kpc, L4 the luminos- 
ity in units of W^Lq, F^fio the flux measured at 60/im in 
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Figure 10. Aio versus K — [12]. Filled and open circles represent 
carbon- and oxygen-rich stars, respectively. The arrowed point, 
the Red Rectangle, is not used in the solution. 



Jansky and Aio the mean wavelength of light emerging from 
the star and its circumstellar dust-shell in units of lO/im. 

Equation 7 assumes all the grains have an emissivity of 
150cm^g~^ at 60/im and a dust to gas ratio of 4.5 x 10~^ 
(Jura 1986). The mean wavelength of light emerging from 
the star and its circumstellar dust-shell, Aio, is given by: 



Al 



\Fxd\ 
r^AdA 



(8) 



where F\ the monochromatic flux density in units of energy 
received per unit area per unit time per unit wavelength 
interval. 

For the stars for which the function Fy could be es- 
timated (section 6.1), Aio was calculated from equation 8 
and is listed in Table 7. The relationship between Aio and 
K — [12] is illustrated in Fig. 10 and can be expressed as: 



Aio = c ± h{K - [12]) ± a{K - [12])^ 



(9) 



where a, b and c are constants. 

Least square flts of equation 9 to the data results in 
the following values (a,b,c) for oxygen-rich (-0.0110,0.284,- 
0.392) and carbon-rich ( 0.0137,0.338,-0.513) stars, respec- 
tively, with an mis error of 0.05 ^m for both chemical types. 
The discrepant carbon star is the Red Rectangle which was 
omitted from the solution. 

Equation 9 with the appropriate parameters is used to 
estimate Aio for the remaining stars in the sample; the re- 
sults are listed in Table 7 

The outflow velocities, wib, were taken to be the expan- 
sion velocities of the shell obtained from C0(1 — > 0) line 
measurements. Appropriate values from the literature are 
listed in Table 7 with references. For the three stars with- 
out measurements an average outflow velocity, for stars with 
Vout < 25kms~^, of 17 ± 4kms~'^ was used. Jones et al. 
(1983, henceforth JHG) have shown there to be a correla- 
tion between outflow velocity and luminosity for a sample of 
OH/IR stars (see flg. 1 in JHG). The outflow velocity versus 
luminosity for stars with periods in our sample is illustrated 
in Fig. 11. It shows no obvious correlation, which should not 
be entirely unexpected, considering the large spread in out- 
flow velocity at a given luminosity (up to 10 km s~^) in flg. 1 
of JHG and the small range in Mm in Fig. 11 (compared to 
that in fig. 1 of JHG). 
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Table 7. Physical and Kinematic data. 
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* - model mass-loss rates taken from Loup ct al. (1993). 

Almost all stars have outflow velocities taken from Loup et al. (1993), with the following exceptions: 
t - outflow velocities taken from Jura & Kleinmann (1989); 
X - average outflow velocity of stars with Vout < 25kms~'^. 
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Log(M„^,i/10-*Mg yr-') 



Figure 11. A comparison of the mass- loss rates determined from 
equation 7 and listed in Table 7 with those given by Loup at al. 
(M^odef) iri the top panel and a plot of outflow velocity versus 
luminosity in the bottom panel (symbols as in Fig. 10). 



The mass-loss rates for 58 of the stars in the sample were 
calculated from equation 7 and are listed in Table 7 together 
with mass-loss rates, Mmodei , taken from the compilation by 
Loup et al. (1993). Where Mmodei was determined from de- 
tailed models of the CO emission. A comparison of these 
values is illustrated in Fig. 11 which shows that they agree 
within a factor of 3. This is reasonable given the uncer- 
tainty in the various assumptions used to derive equation 7. 
17411 — 3154 does not appear in Fig. 11, it has an extremely 
high mass-loss rate, M > lO~*M02/r~^, from equation 7, 
much higher than the value, 0.87 x 10~^MoJ/r~^, given in 
Loup et al. (1993). Kastner (1992) applies a revised model 
of circumstcllar CO emission, and suggests that it is losing 
mass at a rate greater than M > W^^ MQyr~^ , consistent 
with the value derived here. 

Fig. 12 shows the relationship between mass- loss rate 
and K — [12] (M correlates with other colours, but none as 
strongly as this). This relation for Miras has been discussed 
by Whitelock et al. (1994) and Le Bertre & Winters (1998). 
A least squares fit can be made to the data in Fig. 12 as- 
suming the form: 




K-[12l 



Figure 12. Log(M) versus K — [12]. Filled and open circles rep- 
resent carbon- and oxygen-rich stars respectively. The curves are 
described by equation 10. 



AL<1.00 
1.00SAL<1.35 
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Figure 13. Mass-loss rate (M) as a function of pulsation period 



log(M) 



(10) 



(7^-[12]+6) 

where a, b and c axe constants. 

This results in values (a,b,c) for carbon- (-2.89, -1.79, - 
4.32) and oxygen-rich ( 21.34, 3.00, -2.40) stars respectively 
with rms uncertainties of 0.14 and 0.16. This expression can 
then be used to estimate mass-loss rate for the RAFGL 1406, 
which appear in brackets in Table 7. 
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Fig. 13 shows that there is a correlation between mass- 
loss rate and the period of pulsation. Some dependence of 
mass-loss rate on period might be introduced by using the 
PL relation to calculate the quantities L4 and rkpc for equa- 
tion 7. If wo assume, for the sake of argument, that these 
sources are all the same luminosity, then use of the PL 
will result in a mass-loss/period dependence of the form 
M DC P" ""^, i.e. a much weaker effect than is actually seen. 
A correlation between period and mass-loss rate was also 
found by Whitolock ot al. (1994). 

The average mass-loss rate for stars in this sample 
is M ~ 1.5 X 10"^Mo2/T-"^ Stars with P > 1000 day 
(i.e. [12] - [25] > 1.7) have M ~ 1.3 x W'^Msyr-'^ 
and stars with P < 1000 day (i.e. [12] - [25] < 1.7) have 
M ~ 0.87 X lO^'^Moj/r"^ It will be seen later that stars 
with P > 1000 day probably have main-sequence masses 
greater that 4Mo, while those with P < 1000 day have 
main-sequence masses between 1 and 2Mq. Thus interme- 
diate mass stars lose mass an order of magnitude faster 
than do low mass-stars. The average mass-loss rates of 
carbon- and oxygen-rich (including S-type) stars are M ~ 
(l.Oil.O) X IQ-^MQyr-^ and M ~ (2.0±1.2) x W-^Meyr'^ 
respectively. This shght difference may be a consequence of 
assuming the same grain emissivity and gas to dust ratio, in 
equation 7, for the two chemical types. Within the limits of 
the assumptions used here, the mass-loss rates of the carbon 
and oxygen-rich stars are plausibly identical. 

One of the stars in this sample, RAFGL 1406 (lY 
Hya), is a known super-lithium-rich carbon star (i.e. it has 
log > —8) and it is estimated that about 2 percent of 
carbon-stars are of this type (Abia et al. 1993a). Further- 
more there is evidence that these stars are currently impor- 
tant contributors to lithium in the Galaxy (Abia et al. 1995 
and 1998). 

Abia et al. (1993b) estimate a lithium production rate 
of < my 4.5 X 10~^^ M^yr^^ per carbon star with 

a return rate into the interstellar medium P7 = 1.9 x 
10~^ Mqpc~^ Gyr~^. The sample of stars they studied have 
a typical mass-loss rate of 2 x 10~^ MQj/r~^, i.e. more than 
a order of magnitude less than the average mass-loss rate of 
carbon-rich DE-AGB stars. Note that Abia et al. (1993b) 
derive a average progenitor mgiss for the carbon stars in 
their sample of 1.2-1.6 M0, which is similar to the average 
progenitor mass estimated for this sample (see later). As- 
suming a supcr-litiiium-rich carbon star to have a typical 
lithium abundance of log ^^^^ = — 7 (see fig. 3 in Abia et al. 
1993a) we estimate a lithium production rate from carbon- 
rich DE-AGB stars of 1.3 x 10"" M^yr'^ per carbon star, 
with the fraction of super-lithium-rich carbon stars in this 
sample taken to be 0.02. The total return rate of lithium 
by these stars into the interstellar medium is then about 
1.0x10"^ M0pc~^ Gj/r~^. This indicates that carbon-rich 
DE-AGB stars may contribute up to a third of the lithium 
in the Galaxy produced from carbon-rich AGB stars. How- 
ever, granting the underlying assumptions and uncertainty 
in the numbers above, further investigation into the problem 
is necessary (see e.g. Abia et al. 1993b). 



8 SCALE-HEIGHT, PROGENITOR MASS AND 
LIFETIME 
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Figure 14. The function /|^| gives the fraction of stars with 
distances from the plane less than \z\. The curve is described by 
equation 12. 

8.1 Scale- height 

Assuming the usually exponential distribution of stellar den- 
sities the number of stars per height interval, p, varies with 
distance from the galactic plane, z, as: 

p = poe ' ^' , (11) 

where po is the density of stars at a = and H is the scale- 
height. For such a population the average height, \z\, is equal 
to the scale- height, H, provided the sample is reasonably 
complete (as JK89 suggest it is). 

The sample of stars discussed here were necessarily se- 
lected from a limited region of the sky by JK98. Taking this 
into account we use only stars with 14° < I < 230° and find 
that H = 235 pc. 

Alternatively, the fraction of stars, f\^\, within a given 
l^l of the galactic plane is given by: 

/|,| = l-c-l^l/^. (12) 

/j^l is plotted in Fig. 14 for each 100 pc interval for the sam- 
ple of stars with 14° < I < 230° . A curve of the form given 
by equation 12, fitted by least squares, is also shown. This 
fits the data very well (rms difference 0.03) and provides a 
scale-height oi H = 236 ± 10 pc, comparable to the estimate 
made above. This value will be used below as representative 
of the scale height of the sample. 

8.2 Progenitor Mass 

The scale height derived above can be used to estimate a 
typical progenitor mass for the AGB population. A com- 
parison with Miller & Scalo (1979 - their table 1) indicates 
an average main-sequence mass M, ~ 1.2Mo for the group 
under discussion, although it is difficult to estimate the un- 
certainty on this value. 

An alternative approach uses the data in table 3 of 
Miller & Scalo (1979) which gives the number of main- 
sequence stars pc~^(log(M))~^, 0,, as a function of log(M). 
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Figure 15. The interdependence of lower and upper mass limits 
for different scale heights. 



A maximum likelihood straight lino fit to these data over the 
range -0.14 < log(M) < 0.92 has a slope a = -4.20 ±0.32, 
a zero-point /3 = 1.20 ± 0.14 and fits the data fairly well. 
Let Miower and Mapper be the lower and upper progenitor 
mass limits respectively of the stars in this sample. Then 
it is straight forward to show that the average mass, M, in 
Mq is given by: 



M = 



^*-^upper 



a + 1 MSpper - Mf„„,, 
To get an idea of the values of Miower and Mu 



(13) 



the 



relation between Mjou 



and the average scale- 



height H was investigated. H, for a given mass range 



[Mu 



H 



, Mupper] is simply given by: 



M„ 



n(M)(lX(M) 



j^-2';dN{M) 



(14) 



where H{M) is the scale-height as a function of mass M 
(table 1 from Miller & Scalo), dN{M) = jj^^dM the 
number of stars pc^^ in the infinitesimally small interval dM 
centred around M. 

Assuming an average scale- height, Ho, then for a given 
upper mass limit Mupper equation 14 can be solved nu- 
merically for Miower, since the functions H(M) and (fti are 
known. This was done for Ho — H = 236 pc. The resulting 



relation between Miower and M„ 



is plotted in Fig. 15, 



from which it is clear that for Mupper > SMq , Miower has a 
constant value of 0.96Mo. 

Now since we expect to find dust envelope (DE) AGB 
main-sequence progenitors up to about SA^q and since the 
value of M, from equation 13, is relatively insensitive to the 

value of Mupper if ^'^^''^ > 3, values of Miower = O.93M0 

and Mupper = 8Mq wiH^fee assumed. It is important to note 
that the value of A/ is essentially proportional to the value 
of Miower and almost independent of Mupper if '^^'''"^ > 3. 
Also when calculating the surface density of progenitors, it 
will be seen that this is very sensitive to the value of Miower- 



Table 8. Initial Masses and Velocity Dispersions for Main- 
Sequence Stars 

Bp (V) AO AG F5 GO Go 

a (kms^l) 11 14 21 22 23 

M (Mq) 3.2 2.1 1.3 1.1 0.9 



Changing the assumed average scale-height to 216 pc or 
256 pc (i.e. effectively changing the assumed scale-height by 
about 2a) gives values of Miower, assuming Mupper = 8Mq, 
of 0.99Mo and O.86M0 respectively. Thus an error in Miower 
of O.O7M0 will be assumed. This then gives an average mass 
from equation 13 of (1.22 ± 0.18)Mq. 

From the discussion above, an average scale-height of 
236 pc implies a progenitor mass range starting from about 
0.96 Mq up to possibly 8Mq. Hence using equation 18 (be- 
low), one can estimate the ratio of the number of stars in 
the mass range 2.OM0 < M, < 8.0 Mq to the number of 
stars in the mass range I.OMq < Mi < 2.0 Mq. This ra- 
tio is only 0.06. The three stars with periods greater than 
1000 days have an average height above the galactic plane of 
20 ± 2 pc. This indicates that the progenitors of these stars 
were probably more mgissive stars, with Mi > 4.OM0 , as can 
be seen from Miller & Scalo's table 1. Interestingly the ratio 
of stars with P > 1000 day to the rest of the sample (which 
are expected to have P < 1000 day from their colours) is 
0.05, comparable to the ratio calculated above. The ratio, 
from equation 18, of the number of stars in the mass range 
O.93M0 < Mi < 1.0 M0 to the number stars in the mass 
range I.OMq < Mi < 2.0 Mq is about 0.4. Thus the major- 
ity of the stars in this sample have main-sequence masses in 
the range I.OMq < Mi < 2.0 Mq. 

Limited discussion in the literature also supports this 
conclusion. Thus Kahane et al. (2000) estimate the main- 
sequence progenitor of 09452 + 1330, one of the stars in our 
sample, to be of low mass {Mi < 2Mq) by comparing the 
predictions from AGB stellar models of solar composition 
and various initial mass with the observed chlorine isotopic 
ratio in its circumstellar shell. Also Nishida et al. (2000) es- 
timate the progenitor mass of three carbon-rich Miras, with 
P > 450 day, in Magellanic Cloud clusters to be 1.5-1.6 
by main-sequence isochrone fitting. 

The velocity dispersion of the stars with P < 1000 pro- 
vides another indication of initial mass. Fifty three sources 
have velocities, relative to the local standard of rest, tab- 
ulated by Loup et al. (1993), they show a dispersion a = 
24 ± 2kms~^. After correcting for galactic rotation, using 
the constants from Feast & Whitelock (1997), this reduces 
to (J = 19 ± 2kms^^. Table 8 lists the velocity disper- 
sions (calculated from table 2 of Delhaye 1965) and initial 
masses (from Allen 1973) for main-sequence stars between 
AO and G5. The velocity dispersion thus indicates that the 
mean sample has Mj < 2 Mq , consistent with numbers given 
above. 



8.3 Surface Density 

Fig. 16 shows the surface distribution of the stars in the 
galactic plane, calculated from the distances derived previ- 
ously. As can be seen, most of the stars (90 percent) are 
within 1.4 kpc of the sun. There are very few in quadrant 
IV due to positional selection eflfects in the original sam- 
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Figure 16. The distribution of stars around the sun in the Galac- 
tic Plane. Open and closed symbols represent oxygen- and carbon- 
rich stars respectively. 



pie. The surface density, at in region i, i — I, II, III, was 
estimated from: 

a. = ^^ starskpc- 

where A^i is the number of stars in a particular region i. 

The average surface density, (tde-agb of DE-AGB 
stars in the solar neighbourhood Wcis then estimated from: 

2(7/ + an + am , . 

(TDE-AGB = (i-b) 



(15) 



where a J ~ ct/v, assuming symmetry about I = 0. 
Using equations 15 and 16 we find aoE-AGB - 
stars kpc~^. 



15 ± 4 



8.4 Lifetime 

Assuming the average lifetime of stars in any given phase 
of evolution is proportional to the number of stars in that 
phase, it follows that: 

aoE-AGB . 



tde-agb 



ctms 



-tms 



(17) 



where tde-agb is the average lifetime of DE-AGB stars, 
Tms and aMS the average lifetime and surface density, re- 
spectively, of the progenitors of DE-AGB stars on the main- 
sequence. 

From the linear relation between log(^i) and log(M) 
it is straight forward to show that aMS in the mass range 
[Ml Mupper], is simply given by: 

in/f" — M" 1 



aln{W) ' 



(18) 



where <tms is in units of stars pc~ . 

It is clear from equation 18 and the value of a that ctms 
is very sensitive to the value of Miower and a lower if 
j^"'''""' > 3.0. Using the upper and lower mass limits as- 
sumed above, gives aMS = 2.22 ± 0.76 stars pc~^. 



Miller & Scalo's table 1 also gives log(TMs) as a function 
of log(M). Tms is simply given by: 



TMS 



J IV. 



■ TMs{M)dN{M) 



jM^pper ^JV(M) 



(19) 



Equation 19 was integrated numerically with the above mass 
limits, giving tms = (5.5±1.5) x 10® years. The errors above 
were estimated by changing the average scale-height (and 
hence Miower) to 216 and 256 respectively. 

Using equation 17 an average lifetime for the stars in 
this sample weis estimated, giving tde-agb = (3.7 ± 1.9) x 
IC* years. This value is in excellent agreement with lifetimes 
of similar stars at the tip of the AGB estimated by various 
other investigators (see Habing 1996 p. 183). 



9 EVOLUTIONARY STATUS AND 
CONCLUSIONS 

The stars under discussion were selected to isolate high 
mass-loss, thick dust-shell AGB stars within the Solar 
Neighbourhood. This study indicates that the selection cri- 
teria employed were fairly successful in this respect. The 
JHKL photometry, the IRAS two-colour diagrams and 
other published data prove that the majority of these sources 
are cool variable stars with moderate to high mass-loss rates. 

As shown above most of these stars will have had a 
main-sequence progenitor with 1 < M < 2 Mq. While the 
average progenitor mass is ~ 1.3M0. Using this average 
value and an empirical "initial/final mass" -relation (Wei- 
demann 1987) an average final white dwarf mass can be 
estimated at ^ 0.56Mq. 

Assuming the core masses of these stars are close to the 
final white dwarf mass, it is found that the thermal interflash 
period (tif) will be, tij IQp years, using the core-mass 
interflash period relation given by Boothroyd & Sackmann 
(1988) (their equation 6 which is valid for solar composition). 
A more recent core-mass interflash period relation given by 
equation 11a of Wagenhuber & Groenewegen (1998) also 
suggests a thermal interflash period oi Uf ~ 10^ years; this 
is provided there is no hot bottom burning and that the 
star has already undergone several helium shell flashes (see 
Wagenhuber & Groenewegen). The first assumptions is rea- 
sonable as stars under 2 M© do not experience hot bottom 
burning. The second is also an essential assumption for us- 
ing the final white dwarf mass for the current core mass and 
obviously the carbon stars, at least, have undergone several 
flashes. 

This interflash period, tij ~ 10 ', is more than a factor of 
two larger than the average lifetime of these DE-AGB stars, 
Tde-agb = (3.7 ± 1.9) x lO'' years, deduced above. Hence 
we can conclude that these DE-AGB stars will experience at 
most one thermal pulse before leaving the TP- AGB phase. 

The average mass-loss rate, Mave for stars in the sam- 
ple was found to be ~ 1.5 x lO~^M0j/r~^. Thus the to- 
tal mass lost, Miost, by an average star in this phase is 
Miost = Mavef DE-AGB ^ 0.5Mq. Wliicli is vcry close the 
average initial-final mass difference for these stars. However, 
this value for Miost is very uncertain due to the uncertainty 
in the mass-loss rates. Nevertheless, this indicates that a sig- 
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nificant fraction of the envelope mass is lost during this high 
mass-loss dust-enshrouded TP- AGB phase. 

The total rate of mass returned into the interstel- 
lar medium, Mism, by stars in this phase is Mjsm = 
MaveUDE-AGB ~ 2.3 X 10"^ Mq kpc~^ yr~^ . Jura & Klein- 
mann (1990) estimate that M supergiants in the solar neigh- 
bourhood return mass to the interstellar medium at only a 
rate between 1 x lO"'^ and 3 x 10~^ Mq kpc~'^ yr~'^ . Thus 
AGB stars must be important contributors of material to 
the interstellar medium in the solar neighbourhood. 
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APPENDIX 

AVERAGE MAGNITUDES AND COLOURS 
Average Magnitudes 

For variable stars, such as those discussed here, it is often 
useful to know the average magnitude at a particular wave- 
length. The average magnitude as discussed in section 3, was 
taken as the mean of the minimum and maximum magni- 
tude measured for the particular star. This average, Am, is 
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compared in Table Al with the average magnitude, M, de- 
rived from the fit given by equation 4.1 for the 19 stars with 
sufficient data. The average absolute difference, Am, be- 
tween these two averages, for each near-infrared band, with 
the sample standard deviation, a a, on this average is also 
listed in Table Al. From these values for Am it can be seen 
that Am represents the average magnitude fairly well. Ob- 
viously in instances where there are insufficient observations 
to cover a complete cycle. Am will be a poor representation 
of the mean. 



Average Colours 

The stars in this sample not only vary in light output, but 
also in colour. In section 5 the average colour of a star 
is taken as the mean of the value calculated at each Ju- 
lian Date, A. This average is then compared to the average 
colour, Af, derived from the Fourier fits given by equation 
4.1. Af is simply given by: 

Af =< M1-M2 > = < Ml > - < M2 >= Ml,av- M2,av{l) 

where Mi and M2 are the magnitudes in the bands making 
up the colour. Mi,a« and M2,av are the average magnitudes 
in these bands obtained from the Fourier fits given by equa- 
tion 4.1. 

Both the averages, for each ncar-infrarod colour, arc 
listed in Table A2 for the 18 stars with sufficient data. Also 
listed in the Table are the average absolute difference, Aa 
and sample standard deviation, cta, on this average. The 
two values are quite comparable. 
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Table 1. Average magnitudes for variable sources with with good phase coverage 



NAME 




J 






K 




L 






M 


Am 


M 


Am 


M 


Am 


M 


Am 


01037 + 1219 


7.84 


8.06 


4.57 


4.71 


2.34 


2.32 


-0.08 - 


-0.09 


02270 - 2619 


4.54 


4.57 


2.77 


2.74 


1.43 


1.42 


-0.04 - 


-0.06 


02351 - 2711 


2.75 


2.92 


1.57 


1.67 


0.93 


0.98 


0.17 


0.19 


03507+ 1115 


2.23 


2.57 


0.47 


0.69 


-0.59 - 


-0.47 


-1.83 - 


-1.79 


08088 - 3243 


8.45 


8.57 


5.87 


5.97 


3.77 


3.86 


1.19 


1.24 


09116 - 2439 






9.40 


9.68 


6.03 


6.04 


2.15 


2.22 


09429 - 2148 


5.50 


5.65 


3.75 


3.69 


2.40 


2.36 


0.72 


0.73 


09452 + 1330 


7.26 


7.28 


3.98 


3.99 


1.15 


1.17 


-2.54 - 


-2.55 


12447 + 0425 


4.89 


4.88 


3.13 


3.08 


1.83 


1.81 


0.28 


0.26 


17049 - 2440 






8.50 


8.57 


5.47 


5.51 


1.88 


1.93 


17360 - 3012 


8.37 


6.23 


5.98 


6.15 


3.99 


4.09 


1.66 


1.64 


17411 — 3154 














3.86 


3.94 




9.22 


9.33 


6.06 


6.13 


o.u / 




U.O£7 


u.oo 


18333 + 0533 


9.93 


10.38 


5.97 


6.20 


3.57 


3.73 


1.22 


1.29 


18348 - 0526 










8.20 


8.26 


2.01 


1.98 


18398 - 0220 


5.99 


5.91 


3.58 


3.51 


1.87 


1.79 


-0.13 - 


-0.19 


18560 - 2954 


3.03 


3.14 


1.60 


1.67 


0.80 


0.82 


-0.17 - 


-0.19 


19093 - 3256 


3.11 


3.14 


1.88 


1.90 


1.21 


1.22 


0.40 


0.39 


20077 - 0625 


6.70 


6.93 


3.92 


4.04 


2.10 


2.14 


0.18 


0.16 


Am 


0.28 


0.11 


0.05 




0.03 






( 


).1() 


().():! 


0.02 




0.01 





Table 2. Average colours for variable sources with with good phase coverage 



NAME J-H H-K J-K K-L 

(IRAS) Af A Af A Af A Af A 



01037 


+ 1219 


3, 


.28 


3, 


,25 


2.23 


2, 


,21 


5, 


,51 


5.46 


2.42 


2, 


.41 


02270 


-2619 


1, 


.77 


1, 


.76 


1.34 


1, 


,33 


3, 


,11 


3.09 


1.48 


1, 


.47 


02351 


- 2711 


1, 


.18 


1, 


,18 


0.65 


0, 


,65 


1, 


,82 


1.83 


0.76 


0, 


.74 


03507 


+ 1115 


1, 


,77 


1, 


,79 


1.06 


1, 


,08 


2, 


,83 


2.88 


1.24 


1, 


,24 


08088 


- 3243 


2, 


,58 


2, 


,58 


2.10 


2, 


,10 


4, 


,68 


4.69 


2.58 


2, 


.59 


09116 


- 2439 






3, 


,72 


3.37 


3, 


,37 






7.46 


3.88 


3 


.88 


09429 


- 2148 


1. 


,74 


1, 


,74 


1.35 


1, 


,35 


3, 


,10 


3.10 


1.68 


1, 


.67 


09452 


+ 1330 


3. 


,27 


3, 


,28 


2.83 


2, 


,83 


6, 


,10 


6.11 


3.69 


3, 


.70 


12447 


+ 0425 


1. 


,76 


1, 


,75 


1.30 


1, 


,29 


3, 


,06 


3.04 


1.55 


1, 


.55 


17049 


- 2440 






3, 


,24 


3.03 


3, 


,04 






6.20 


3.60 


3 


.62 


17360 


-3012 


2. 


,39 


2, 


,99 


1.99 


1, 


,99 


4, 


,38 


4.98 


2.33 


2, 


.34 


18194 


-2708 


3. 


,16 


3, 


,16 


2.40 


2, 


,39 


5, 


,55 


5.54 


2.78 


2, 


.77 


18333 


+ 0533 


3. 


,96 


3, 


,89 


2.40 


2, 


,34 


6, 


,35 


6.24 


2.36 


2, 


.32 


18348 


- 0526 






















6.20 


6, 


,27 


18398 


- 0220 


2. 


,40 


2, 


,37 


1.71 


1, 


,68 


4, 


,12 


4.05 


2.00 


1, 


.98 


18560 


- 2954 


1. 


,44 


1, 


,46 


0.79 


0, 


,82 


2, 


,23 


2.28 


0.98 


1, 


.00 


19093 


-3256 


1. 


,23 


1, 


,22 


0.67 


0, 


,65 


1, 


,90 


1.87 


0.81 


0, 


.77 


20077 


-0625 


2. 


,79 


2, 


,76 


1.82 


1, 


,79 


4, 


,60 


4.55 


1.92 


1, 


.91 



Aa 0.06 0.01 0.07 0.02 

0.04 0.01 0.04 0.01 
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